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Abstract 



Most of the exact solutions of quantum one-dimensional Hamiltonians are 
obtained thanks to the success of the Bethe ansatz on its several formula- 
tions. According to this ansatz the amplitudes of the eigenfunctions of the 
Hamiltonian are given by a sum of permutations of appropriate plane waves. 
In this paper, alternatively, we present a matrix product ansatz that asserts 
that those amplitudes are given in terms of a matrix product. The eigenvalue 
equation for the Hamiltonian define the algebraic properties of the matrices 
defining the amplitudes. The existence of a consistent algebra imply the ex- 
act integrability of the model. The matrix product ansatz we propose allow 
an unified and simple formulation of several exact integrable Hamiltonians. 
In order to introduce and illustrate this ansatz we present the exact solu- 
tions of several quantum chains with one and two global conservation laws 
and periodic boundaries such as the XXZ chain, spin-1 Fateev-Zamolodchikov 
model, Izergin-Korepin model, Sutherland model, t-J model, Hubbard model, 
etc. Formulation of the matrix product ansatz for quantum chains with open 
ends is also possible. As an illustration we present the exact solution of an 
extended XXZ chain with z-magnetic fields at the surface and arbitrary hard- 



core exclusion among the spins. 



I. INTRODUCTION 

The Bethe ansatz [1] and its generalizations merged along the years as a quite effi- 
cient and powerful tool for the exact solution of the eigenspectrum of a great variety of 
one-dimensional quantum chains and two-dimensional transfer matrices (see e.g. [2]- [5] for 
reviews). According to this ansatz the amplitudes of the wavefunctions are expressed by 
a nonlinear combination of properly defined plane waves. Although the question of com- 
pleteness of the Bethe solutions is in general an open and difficult problem [6,7], a quantum 
Hamiltonian is considered exactly integrable if, in the thermodynamic limit, an infinite 
number of its eigenfunctions are given by the Bethe ansatz. 

On the other hand along the last two decades [8]- [11], several models were discovered 
having a ground state wavefunction exactly known and obtained though an ansatz known as 
the matrix product ansatz. Differently from the Bethe ansatz, this approach only gives the 
ground state eigenfunction, whose amplitudes are expressed in terms of a product of matrices, 
or more generally in terms of a product of generators of quadratic algebras [12]. In a distinct 
context a matrix product ansatz has also been applied quite successfully to the exact solution 
of the stationary distribution of probabilities of some one-dimensional stochastic models [13]. 
The similarity between the master equation describing the time fluctuations of these models 
and the Schrodinger equation in Euclidean times enables to identify an associated "quantum" 
Hamiltonian to these stochastic models. The simplest example is the problem of asymmetric 
diffusion of hard-core particles on the one-dimensional lattice (see [14]- [17] for reviews). 
The time evolution operator governing the time fluctuations of this last model coincides 
with the exact integrable anisotropic Heisenberg chain, or XXZ chain, on its ferromagnetic 
regime and appropriate boundary fields [17]. The related quantum chains in general are 
not exact integrable but have their ground state eigenfunctions given in terms of a matrix 
product ansatz. The matrix product ansatz, although providing only stationary properties of 
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some stochastic system produced interesting results in a quite variety of problem including 
interface growth [18], boundary induced phase transitions [13], [19]- [21] the dynamics of 
shocks [22] or traffic flow [23]. 

More recently an interesting development of the matrix product ansatz, named dynam- 
ical matrix product ansatz [24], merged in the area of stochastic one-dimensional systems. 
Models satisfying this ansatz, distinctly from the previous formulated matrix product ansatz, 
have their probability distributions, at arbitrary times, given in terms of a product of ma- 
trices, which are now time dependent. The dynamical matrix product ansatz was shown 
originally to be valid in the problem of asymmetric diffusion of particles on the lattice [24,25]. 
More recently [26,27] the validity of this ansatz was also confirmed in the exact integrable 
manifold of the asymmetric diffusion of particles with two kinds of particles. The validity 
of this ansatz for such integrable system induces us to expect that all the quantum chains, 
related or not to the stochastic systems, that are solvable though the Bethe ansatz may also 
be solvable by an appropriate matrix product ansatz. We expect that all the components 
of an arbitrary eigenfunction of an exact integrable quantum chain, that are normally given 
in terms of a combination of plane waves, can also be expressed in terms of a product of 
matrices satisfying algebraic properties that ensure the exact integrability of the model. 

In this paper we are going to show the validity of this conjecture for a huge family 
of exact integrable quantum Hamiltonians, by showing how to formulate their solutions 
in terms of a matrix product ansatz. A brief summary of some of our results has been 
announced in [28]. We are going to present explicit examples of a matrix product ansatz 
formulation for models having one or two global conservations. Examples of models with 
one global conservation law (U(l) symmetric) include the XXZ quantum chain [29], the 
spin-1 Fateev-Zamolodchikov model [30], the Izergin-Korepin [31] and the solvable spin-1 
model considered in [32]. Among the models with two global conservation laws (U(1)®U(1) 
symmetric) we present the matrix product ansatz for the models where, in its formulation 
in terms of particles, no double occupancy of sites is allowed, like the anisotropic spin-1 
Sutherland model [33], the spin-1 Perk-Schultz models [34], the fermionic supersymmetric 
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t-J model [35], and also for the models with double occupancy like the spin-3/2 Perk-Schultz 
model [34], the Essler-Korepin-Schoutens model [36], the Hubbard [37] model and the two- 
parameter integrable model introduced in [38]. 

The ansatz we produce enable us to formulate in a simple and unified way generalized 
integrable models where the quantum spins have now hard-core interactions that exclude the 
occupation of two spins at neighbouring sites. Although most of our results will be presented 
for periodic chains they can also be extended for non periodic but integrable chains. As an 
example of a non periodic chain we present the solution, in terms of a matrix product ansatz, 
for the exact integrable XXZ chain with boundary z-magnetic fields [39]. 

The layout of the paper is as follows. In sections 2 and 3 we consider models with a 
single global conservation law. In section 2 we introduce and solve through our matrix 
product ansatz a generalization of the standard XXZ chain, where the up spins have a 
hard core interaction of a given, but arbitrary, range in terms of units of lattice spacing. 
In section 3 models with spin 1 are considered. The matrix product ansatz is formulated 
for the spin-1 Fateev-Zamolodchikov model, the Izergin-Korepin model, as well as for the 
spin-1 model introduced in [32]. All these solutions are obtained through a single matrix 
product ansatz, whose matrices satisfy for each model distinct algebraic relations. In section 
4 and 5 we consider models with U(1)®U(1) symmetry, thus exhibiting two conserved global 
quantities. In section 4 we consider the models of spin-1 in this class. Those models are the 
anisotropic Sutherland model, or the SU(3) Perk-Schultz model and the supersymmetric t-J 
model. In section 4 we present a general solution that includes the spin-3/2 Perk-Schultz 
model, the Essler-Korepin-Schoutens model, the Hubbard model as well the general exact 
integrable two-parameter model presented in [38]. Again the solution of all these models 
are presented through an unified matrix product ansatz. All the solutions presented in the 
previous sections are derived for quantum chain defined on lattices with toroidal boundary 
conditions, in section 6 we show how to extend our solutions for the case of the XXZ chain 
in an open lattice with a z-magnetic field at the end points of the lattice. Finally in section 
7 we close our paper with some conclusions and comments. 

4 



II. GENERALIZED XXZ CHAINS 



As a first application of our matrix product ansatz we are going to present, in this section, 
the exact solution of the anisotropic Heisenberg model or XXZ chain. In order to show the 
powerfulness of our ansatz, instead of solving the standard XXZ chain we are going to solve a 
generalization of this quantum Hamiltonian where arbitrary excluded volume are considered 
among the up spins (<r z -basis). We consider a generalized XXZ chain where any two up spins, 
due to hard-core interactions, are not allowed to occupy lattice sites at distances smaller 
than s (s = 1, 2, . . .), in units of lattice spacing. Unlike the down spins, that have only on 
site hard-core exclusions, the up spins behave as they would have an effective size s. The 
generalized XXZ Hamiltonian we consider, in a lattice with L sites, is given by 



where a x , a y , a z are spin-1/2 Pauli matrices, A is the anisotropy and the projector V s 
projects out from the associated Hilbert space the configurations where any two up spins 
are at distance smaller than s. In the particular case where s = 1 the projector V s is the 
unity operator and we recover the standard XXZ chain. The Hamiltonian (1), in terms of 
raising and lowering spin-1/2 operators a ± = (a x ±ia v )/2 is a particular case of the general 
Hamiltonian 



where e+ = e_ = 1 and a harmless constant (+LA/2) was added. Without any additional 
difficulty we are going to consider the solution of (2), for general values of e+, e_ and A. If 
we interpret the up spins as particles and the down spins as vacant sites the Hamiltonian (2) 
with the choice e + + e_ = 1 = — 2 A coincides, apart from a harmless constant, with the time 
evolution operator of the asymmetric diffusion problem (asymmetry e + /e_) of particles with 
size s on the lattice [40]. We are going to consider (2) with the periodic boundary condition 




i=l 



(1) 





<J L+ i = o x , 



_z _z 
°L+1 = °\- 



(3) 
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In section 6 we are going to consider the case of open boundary conditions. 

The Hamiltonian (2) with the boundary condition (3) has a XJ(l)<g>Z L symmetry due 
to its commutation with the total spin operator S z = Yh=\ ° z an d the spatial translation 
operator T = e lP on the lattice. Consequently the Hilbert space associated to (2) can be 
separated into block disjoint sectors labelled by the number n of up spins (n — 0, 1, . . . , L) 
and the momentum eigenvalues P (P — ^Z, I — 0, 1, . . . , L — 1). 

The ansatz we propose asserts that any eigenfunction \^ n ,p) of (2) in the sector with 
n up spins and momentum P, will have its components given by the trace of the matrix 
product 

l*n,p) = J2 f(xi,...,X n )\x 1 ,...,X n ), 

f(x u ...,x n )=Tr [E Xl - 1 A^E Xa ~ Xl - 1 A^ ■ ■ ■ E^'^-^A^ E L ~ Xn n P ] . (4) 

The configuration in cr^-basis where the up spins are located at (x 1 , . . . ,x n ) are denoted 
by the ket | ) and the symbol (*) in the sum denotes the restriction to the sets 

satisfying the hard-core exclusion due to the effective size s of the up spins, i.e., 

x i+i > Xi + s, i — 1, . . . , n — 1, Xi > 1, s < x n — xi < L — s. (5) 

Differently from the standard Bethe ansatz where f(x\, . . . , x n ) is given by the combination 
of plane waves now it is given by the trace of a product of matrices. Actually E, A^ 
and flp are abstract objects with appropriate commutation relations. These operators in 
order to provide well defined amplitudes for the eigenfunctions should have an associative 
product, and as usual in the literature, we call them simply as matrices. The matrices E 
and A^ are associated to the down and up spins representing the ket configuration. The 
superscript s is just to remember that the up spins have an effective size s. The matrix flp 
in (4) is introduced in order to fix the momentum P of the eigenfunction |^ n; p). This is 
accomplished by imposing the algebraic relations 1 



L In the general case of twisted boundary conditions cr^+i = e±i * CJ i : ) = °i the ansatz also 



6 



(6) 



since from (4), eigenf unctions of momentum P should have the ratio of amplitudes 

. . . , x n ) 



-iPm 



(m = 0,l,...,L-l). 



(7) 



f{xi + m, . . . ,x n + m) 

The matrix product ansatz for \^ n ,p) will work if we are able to obtain consistent algebraic 
relations among E and that solve the eigenvalue equation 



H 8 \V n , P ) = S n \^n,p). 



(8) 



As is customary, before considering the case of general values of n let us consider the 
eigensectors with n — 1 and n = 2 up spins. 

n = 1. For one up spin the eigenvalue equation (8) give us for the amplitudes with a 
single spin at X\ — 1, . . . , L 



eiTr 



E xi - 1 A( s) E L - xi n, 



-e+Tr 
+ATr 



E xi - 2 A is) E L - xi+1 n P 

E xl - l A {s) E L - Xl VLp 



-e_Tr 



E x lA {s) E L- Xl -l n 



(9) 



The cyclic property of the trace and relation (6) fix the eigenenergies 

e 1 = - (e + e~ iP + e_e lP - A) , 



(10) 



where P = ™ (I — 0, 1, . . . , L — 1) is the momentum of the state. An alternative solution 
of (9) that will be simple to generalize for arbitrary values of n is obtained by replacing 



A^ = A { k s) E 2 - s , 



(11) 



where now A^ is a spectral parameter dependent matrix obeying the following commutation 
relation with the matrix E: 



EA ( k s) = e lk A ( k s) E. 



(12) 



works with an appropriate generalization of (6). 
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Inserting (11) in (9) and using (12) we obtain 



£l = £ (k) = - [e + e~ ik + e_e jfc - A) , 



(13) 



and comparing (10) with (13) we fix the spectral parameter k as the momentum of the 
1-particle eigenfunction \^i,p), i. e., k = P = ™ (I — 0, 1, . . . , L - 1). 

n = 2. The eigenvalue equation (8) gives for the amplitudes of \^ 2 ,p) two types of 
relations depending if the two up spins located at x\ and x 2 are at the closest position or 
not. The amplitudes corresponding to the kets \x±, x 2 ) where x 2 > X\ + s give us the relation 



s 2 Tr r^i- 1 ^ E X "- X1 - 1 A^ E L ~ X2 Q P 
-e_Tr \e Xi A (s) E X2 - xl ~ 2 A is) E L ' X2 Q P 
-e_Tr \E xl - 1 A ia) E X2 - xl A w E L - X2 - 1 n P 



-e+Tr 



E xi-2 A (s) E X2-x lA (s) E L-x 2 Q 



e+Tr 



E Xl ~ 1 A^ E X2 ~ X1 ~ 2 A^ E L ~ X2+1 Q) 



+ 2ATr 



E Xl ~ l A^ E X2 ~ X1 - 1 j\( s ) E L ~ X2 Vti 



(14) 



A possible and convenient way to solve this equation is obtained by a generalization of (11) 
and (12). We identify A^ as composed by two spectral parameter dependent new matrices 



4? and 4?, i- e., 



A^ = J2A^E 2 - S , 
i=i 

that satisfy the commutation relations 

EA% = e*>A%E, (4J) 2 = 0, J = 1,2. 
Inserting (15) in (14) and using (16) we obtain 

e 2 = e(fci) +e(k 2 ), 



(15) 



(16) 



(17) 



where e(k) is given by (13). The relation (6) give us the commutation of the new matrices 
4f w ith f2p, i. e., 



A^Q P = e^Q P Ajg j = 1,2. 



(18) 



Comparing as in (7) the amplitudes of the configurations \xi, x 2 ) and \xi + m, x 2 + m) and 
exploring the cyclic property of the trace we obtain 
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P = h + fc 2 . 



(19) 



Up to now we have at our disposal, for solving the eigenvalue equation, the commutation 
relation of and as well the spectral parameters k\ and k 2 , that might be in general 
complex numbers. These commutations relations are going to be fixed by the application of 
the eigenvalue equation (8) to the components of the configurations \xi,x 2 ) where the up 
spins are at the "matching" condition x 2 = X\ + s: 



e 2 Tr 



exi-i^s) E s-i A (s) E L -^- s n, 



-e_Tr 



E xi-l A (s) E s A (s) E L-xi-a-lQ 



= -e+Tr 
+ ATr 



- EX1 -2 A (s) E s A (s) E L-xi-s^ p 

E X!-1 A (S) E S-1 A (S) E L-XlSQ p 



(20) 



Using (15), (17) and the commutation relations (16) in the last expression one obtains 

2 



£ [e_ + e+e-^^ - Ae 

3,1=1 



A ( S ) A ( S ) 
kj k[ 



o. 



(21) 



This imply that the matrices A k : (j = 1, 2) should obey 

AilA'hSik^k^A^, (Z#), (A^f = 0, (j,/ = 1,2), (22) 



where 



S(k j ,k l ) = 



(23) 



This last relation in the context of (1 + 1) -dimensional field theory is know as the Zamolod- 
chikov algebra [41,42]. The complex spectral parameters kj (j = 1,2) that are still free up 
to now are fixed by imposing that the ratio of two components f(xi,X2)/f(x' 1 ,x / 2 ) should 
be uniquely related. The cyclic property of the trace with the algebraic relations (16), (17) 
and (22) give us 



Tr 



A%A%E L - 28+2 n P 



= e -^s + 2)k jTl \ A (°) E L-2s + 2 A (s) n 



= e - ifej L e i2fc,( S -l) e -iP( S -l)£^ 



A^A^E L - 2s+2 n P 



, (24) 



or equivalently, since P — k\ + k 2 , 



, J = 1,2 (j^l). 



(25) 



The energies £2 and momentum P of |^2,p) are obtained by inserting the solutions of (25) 
into (17) and (19) respectively. 

General n. The previous calculation can be easily extended for arbitrary values of the 
number n of up spins. The eigenvalue equation (8) when applied to the amplitudes of \^ n ,p) 
corresponding to the configurations where all the n spins are at distances larger than the 
excluded volume s, gives a generalization of (14): 

e n Tr [• • • e^-^- 1 - 1 A {s) e^-^A^ ■ ■ ■ A (s) E L - Xn Vt P = 

n 

— ^{ e+ Tr ■ • ■ E Xi ~ Xi - 1-2 A^ E Xi+1 ~ Xl A^ ■ ■ ■ A^ E L ~ Xn Q, P 

i=l 

+e_Tr [• • .^i-^-i-^W^+i-^-2^(») . . . A(s)E L - Xn+1 Slp 



— ATr 



}• 



(26) 



The solution of this last equation 2 is obtained by identifying the A^ matrix as a combination 
of n spectral parameter dependent matrices {A^ • j = 1, . . . , n}, i. e., 



1 1 



(27) 



3=1 



with the commutation relations with the matrices E and Q P given by 

EA^ = e^A^E, A$Q P = eW-)n P A$ (j = 1, . . . , n). 



(28) 



The energy and momentum are obtained by inserting (27) into (26) and using (28), and are 
given by 



3=1 3=1 



(29) 



2 The most general relation = Y7j=i E a A^E 13 could be used. However (27) is more convenient 
since otherwise the structure constants S(k{, kj) in (23) will depend on the size s. 
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respectively. The eigenvalue equation (8) gives for the configuration where a pair of up spins 
are at the "colliding" positions x i+1 — Xi + s a relation that coincides with (21)-(23) but 
with j,l — 1, . . . , n. The relation coming from the configuration where three particles are at 
the colliding position Xj+i = X{ + s, Xi + 2 = Xi+i + s is given by 

n n 

E + e + e-*to + * - Ae-^) E 4^ 

j,l=l t=l 
n n 

+ J2 e~ ikj A { £ E ei(kl+kt) { e - + e+e~ l(fc;+fet) - Ae"* 1 ) A% 4? = °, ( 30 ) 

being a consequence of the relation (21). Similarly the amplitudes with arbitrary number 
of particles at colliding position will be automatically satisfied if the matrices E, (j = 
1, . . . , n) and f2p obey the algebraic relations (22) and (28). The associativity of the algebra 
provide a well defined value for any product of matrices and follows from the fact that the 
structure constants S(ki, k 3 ) in (22)-(23) are c-numbers with the property S(ki, kj)S(kj, fcj) = 
1 = l,...,n). 

The cyclic invariance of the trace in (4) will fix the complex spectral parameters k 3 - 
(j = 1, . . . , n), providing a well defined value for the components of l^n.p), i- e., 

Tr [4 S) • • • 4 S) • • • Ai s) E L ~ n ^n P 

n 

= e ^[W.-i)] e «-(i-) JJ S (k 3 , k t )Tr [A%A$ ■ ■ • ^ • • • A^E L ~^n P 

l=j+l 

= e^IW-DIe^d-) JJ S(kj, h)Tr [A$ ■ ■ ■ 4? • • • A^E^'^Qp] , j = 1, . . . , n, 

(31) 

or equivalently 



s-1 

e ^ = (_l)«Jj(£2M c+ " rc "". r ." j = l,...,n. (32) 

The acceptable set {A^-; j — 1, . . . ,n} of spectral parameters defining the eigenvectors |^ n ,p) 
are the solution of (32) where fcj 7^ fcj = 1, . . . ,n). Since (4 s "*) 2 = 0; solutions of (32) 
with coinciding roots give us null states. 

The equation (32) coincides with the Bethe ansatz equation derived for the model (2) 
through the standard Bethe ansatz [40]. The choice e + = e_ = 1 with s = 1 gives the Bethe 
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ansatz equation of the standard XXZ chain (1) [29]. Moreover, using (27), an arbitrary 
component f(xi,..., x n ) of the eigenfunction |^ n ,p), given by the ansatz (4), can be written 

as 



f(x±, . . . , X n ) 

n n n 

EE-ETr 

il=lZ2 = l »n=l 



The commutation relations (28) allow us to write 



(33) 



f{X\ ■>•••■> Xn) 



n n 



^ ... ^ e »[fcn(^i-l)+-+^ n (^n-l)] Tr 

«1 = 1 «n = l 

Let us define the new matrices 



4° £ 1_s 4° • • • £ 1_s 4 s) e l q,p 



(34) 



4 ) = 4 ) ^ S (J = l,-,n). 



(35) 



It is simple to verify, from (22), that they satisfy, for j, I — 1, . . . , n, 



4: ) 4: ) = ^,^)4l ) 4 ) ' ovo, (4?) 2 = °> 



(36) 



where 



s-1 



5(%,A; Z ) = S(%,^) ( 



In terms of these new matrices, and exploring the fact that (4^) 2 = we can write 



(37) 



f(xi, . . . , x n ) — ^2 e i ' fcpi ( xl-1 ) + "' +fcp "( x "~ 1 ^Tr 

pi,—,Pn 



(38) 



where the sum is over the permutations {pi,P2, ■ ■ ■ ,Pn} of the non repeated integers 
(1, . . . ,n). This last result show us that the amplitudes derived using the proposed matrix 
product ansatz (4) is given by a combination of plane waves with complex wave numbers 
{ki}, and reproduces the results previously obtained for the Hamiltonian (2) [40] through 
the coordinate Bethe ansatz. 
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III. MODELS OF SPIN 1 WITH ONE CONSERVATION LAW 



We present in this section the appropriate matrix ansatz solution for models of spin 1 
with a single global conservation law, like the XXZ chain of last section. Integrable models 
on this category are the Fateev-Zamolodchikov quantum chain [30], the Izergin-Korepin 
model [31] and the spin-1 model introduced in [32]. 

We start with a general spin-1 model with U(l) symmetry and nearest neighbor inter- 
action. Instead of writing this general model in terms of spin-1 Pauli matrices it is more 
convenient to write it in terms of the 3x3 Weyl matrices E l ' m (l,m = 0,1,2), with i,j 
elements (E l ' m )ij = 5ij5 m j. At each lattice site we may have zero particles (s z = — 1), 
one particle (s z = 0) or two particles (s z = 1). The general Hamiltonian we consider, that 
conserves the number of particles, is given by 

L 

rjU(l) _ V- , rrOO 
n s=l — ~ 2^ n i,i+l + LlL 0' 
i=l 

E KLK'^i, (< = i,...,l) (39) 

k,l,m,n=0 

where periodic boundary conditions are imposed and in order to conserve the total number 
of particles T^ n = if k + 1 ^ m + n. 

The total number of particles n (0,1, ... , 2L) and the momentum P (^, 1 = 0,..., L—l) 
are good quantum numbers for the general Hamiltonian (39). The eigenfunctions on these 
eigensectors are given by 

* 

l*n,p) = E f(x 1 ,...,X n )\x 1 ,...,X n ) (40) 

where (x\, . . . , x n ) are the coordinates of the particles and the symbol (*) in the sum means 
restriction to the sets with coordinates {x i+ \ > x iy x i+2 > 

In order to formulate an appropriate matrix product ansatz for (39) we associate, as 
in the previous section, the matrices E and A to the empty sites (s z = —1) and the sites 
occupied by a single particle (s z = 0), respectively. The sites with double occupancy of 
particles (s z = 1) are associated with the matrix BE~ X B. Certainly the Hamiltonian (39) 
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is not exact integrable for arbitrary values of {T k rn l n }. The ansatz we propose states that in 
the exact integrable manifold of (39) any eigenfunction |^ n ,p), in the sector with n particles 
and momentum P will have components given in terms of traces of a matrix product. The 
amplitudes corresponding to the configurations where there is no double occupancy are given 
by 



f{x u ...,x n )=Tr \E X1 - 1 AE X *- X1 - 1 A ■ ■ ■ E Xn ~ Xn ~ 1 ~ 1 AE L ~ Xn Qp 



while if there exists a double occupancy at Xj+i = X{ we have 



(41) 



f{%li ■ ■ ■ i %ii ■ ■ ■ i %n) 



Tr 



E Xl ~ l A ■ ■ ■ fiXi-vi-i- 1 b E~ l B E Xi+2 ~ Xi ~ 1 ■ ■ ■ E 



Xn-Xn-l-l A TPL-X 



AE^~ Xn fl P 



(42) 



As in previous section the matrix flp and the trace operation are introduced in order to fix 
the momentum P = ^7, (I — 0, 1, . . . , L — 1) of \^ n ,p), this is accomplished (see (7)) by 
imposing 



EVtp = e- iP tt P E, AVtp = e'^ttpA, BVt P = e'^QpB. 



iP ( 



-iP ( 



(43) 



The algebraic relations among the E, A and B operators will be obtained from the require- 
ment that (39) with the ansatz (40)-(42) are solutions of the eigenvalue equation 



H^l ) \V nt p)=e n \V nt p). 



(44) 



Let us consider initially the cases with a small number of particles, 
n = 1. The eigenvalue equation gives for the amplitudes with a particle at position 
Xi — 1, . . . , N: 



£iTr 



-r 1 °Tr 
1 o 1 11 



E Xl ~ 2 AE L - Xl+1 n, 



- T° x Tr 
1 1 o 11 



E X1 AE L - xl ~ l Vt } 



+ (2rn-rl°-r°i)Tr 



E Xl ~ l AE L - Xi n, 



(45) 



As in previous section a convenient solution is obtained by introducing the spectral param- 
eter matrix Aj,, 
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A = A k E with EA k = e ik A k E. 



(46) 



Inserting (46) into (44) we obtain 



ei = e{k) = - (i* ?e- fe + T? J e * - 21* g + I* ° + C) , 



P = k 



2nl 



(Z = 0,...,L-1). 



(47) 



n = 2. In this case the eigenvalue equation (43) produces relations that are distinct if 
x 2 > X\ + 1, x 2 — Xi + 1 or x 2 — X\. The relations coming from the amplitudes where the 
particles are not at the "colliding" positions x 2 = X\ oi x 2 = X\ + 1 are just a straightforward 
generalization of (45), for two particles, whose solution is obtained by identifying the matrix 
A as composed by two spectral parameter dependent matrices A kl , A k2 , i. e., 



A = ^2A k .E 



(48) 



satisfying the commutations relations 



EA k .=e ik 'A k .E, (A k .) 2 = j = l,..., 



n. 



(49) 



The energy and momentum are given by 



3=1 3=1 



(50) 



with e(k) given by (47). The eigenvalue equation (44) gives for the components f(x 1 , X\ + 1) 
the equation 



£ 2 Tr 

-r° x Tr 
1 1 o 11 

-r??Tr 



E Xl - 1 AAE L - xl ~ 1 Qp 

EX1 -l AEAE L- X1 -2 

E Xl BE- x BE L ~ Xl - x Sl 



-I^Tr [E Xl - 2 AEAE L - Xl - 1 Sl P 

- r?°Tr [E X1 - 1 BE- 1 BE L - X1 SI P 
j- fir 00 r 10 r 01 r n Wv 

+ poo 1 1 1 1 1 1 1) 11 



E Xl ~ x AAE L ~ Xl ~ x Vti 



(51) 



while the components f(xi,xi) gives 
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e 2 Tr 



E xl ~ 1 BE~ 1 BE L ~ xl ^li 



— F 1 lr TY 
1 02 11 



— r 2 °ty 

i 02 il 



E xl ~ 2 AAE L ~ xi n, 



— F 1 lr TY 

1 2 11 



E xl - 1 AAE L - xl - 1 n P 



E Xl ~ 2 BE- l BE L ~ Xl+l Sl } 



1 2 1J: 



E xl BE~ 1 BE L ~ xl ~ 1 n t 



+ (2r° ° - r 2 ° - r° 2 ) Tr [e^be^be^vlp 



(52) 



Both equations (51), (52) are solved by using (48) and by writing the matrix B also as a 
combination of two spectral parameter dependent matrices B kl , B k2 , i. e., 



2 

B = Y,B k] E with EB kj =e ik *B kj E, (B k .) =0 (j = l,...,n). 



(53) 



Notice that the matrices A kj , B k . have the same set {k\, k 2 } of undetermined complex 
spectral parameters. Inserting (48)-(50) and (53) in (51) and (52) we obtain 

n n 

£ ^(%, fc,)Tr [e x * A k .A kl E L -*^ P ] = ]T A^Tr [E^B k .B kl E L -^Q P ] (54) 

and 



3,1=1 



E xl B kj B kl E L - xl tt 



= £ C 2 (k j ,k l )e ik 'Tr 
3,1=1 



E x 'A kj A kl E L ~ Xi n P (55) 



where 



Ci(fcj, fcj) = r 2 ? + r? 2 e i(fcj ' +fc ' ) 

Co(%, fc,) = Tl\{e ih + e ifc + r? J(e ifc ' + e^)e l{fci+fcj) 
, / 9 -pi o , 9 r o i _ r 2 _ r° 2 _ or °^ J( k 3+ k i) _ r 2 _ r° 2J2(kj+h) 

' y 1 1 i 20 i 02 ZL 00) e i 02 1 2 e 

c 2 (fc J -,fc,) = rSJ + r^ e < ^ + *«) 



The relations (54) and (55) imply 



7V(%, /cz)^A + fcj)A fcj i4 fc . = C!(%, fc,) (B k .B kl + B fcl B fcj ) 



(56) 



(57) 



C (kj, h) (B kj B h + B fc ,fl fcj .) = C 2 (kj, h) (e ikl A kj A h + e ik *A kl A kj ) , 



(58) 
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with k, I = 1, . . . , n. Multiplying (58) by the symmetric function Ci(kj, k\) and using (57), 
the last equation gives the relation 



A kj A kl = S{kj, ki)A kl A kj (j I), 

c/k k \ C (kj, k)N(k h kj) - Ci(kj, kjCzjkj, h)e^ 

The spectral parameters k\ and k 2 are going to be fixed by the cyclic property of the traces 
defining the amplitudes (41) and (42). Since the B matrices only appear in the combination 

BE~ 1 B= (jr B kj B k )j E. (60) 

The use of (57) allow us to replace in any amplitude the {B k } matrices by the {A k } matri- 
ces. Then an arbitrary amplitude should be proportional to Tr A kl A kj E L Qp . The cyclic 
property of the trace and the commutation relations (43), (49) and (59) give us 

e^^Sik^h), j = 1,2 (jVO (61) 

with S(kj, ki) given by (59). It is interesting to mention that we obtained for n = 2 particles 
the integrability for arbitrary values of the coupling T^ n of the Hamiltonian (39). Certainly 
this will be not the case for n > 2 particles. 

n = 3. In this case we have several distinct types of relations for the amplitudes 
f(xi, x 2 , x 3 ) in (41) and (42). The eigenvalue equation when applied to the amplitudes 
where all the particle are not at "colliding" positions, i. e., x 3 > x 2 + 1 > x\ + 2 give us 
a straightforward generalization of (45), for three particles, whose solution is obtained by 
identifying, as in (48), the matrix A as composed by three spectral parameter dependent 
matrices A kl , A k2 , A ks . The energy and momentum are given by (50) with n = 3. The 
components f(xi, x 2 = Xi, x 3 ), f(x 1 , x 2 = x 1 + 1, x 3 ) with x 3 > x 2 + 1 and f(xi, x 2 , x 3 = x 2 ), 
f(xi,x 2 ,x 3 = x 2 + 1) with x\ < x 2 — 1 give us generalizations of (57) and (58) for n = 3. We 
have new relations when the eigenvalue equation is applied to the amplitudes corresponding 
to the 3 particles at the colliding positions. This happens for f(x,x,x + 1) 
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e 3 Ti 



E X - 1 BE- 1 BAE L - X - 1 Q, 



— r 1 lr r> 
1 2 11 



— F 2 °Tr 

i 02 il 



E x ~ 2 AAAE L - x - l Vtp 



E x - l BE- l BEAE L - x - 2 tt P 



E X ~ 2 BE~ 1 BEAE 



L-x-l 



E x ~ l ABE- x BE L - x - x Slp 



i /o-nOO F 2 r 2 1 plO\rp 
00 1 2 1 2 1 1 1 0,1 11 

and for amplitudes f(x,x+ l,x+ 1) 



E X ~ 1 BE~ 1 BAE L ~ X ~ 1 Q P 



(62) 



e 3 Tr 



E x ~~ l ABE~ l BE L ~ x ~ 1 



1 2 11 



_r° 2 Tr 

1 2 11 



1 1 - 11 



E x ~ 2 AEBE- l BE L - x - l Sl } 



E x ~ l AEBE- l BE L ~ x - 2 VL } 



E x ~ l BE~ l BAE 



L-x-l 



/o-nOO _ F 2 _ r l 2 _ r lA T 
00 i 20 i 12 L 01J 11 



E X ~ 1 ABE~ L BE 



1 D ETii— X— 1 



(63) 



Inserting (48), (50) and (53) with n = 3 in (62) and (63) we obtain the new algebraic 
relations relating three spectral parameter matrices. Using (57) and (58) to simplify those 
expressions we obtain 

]T (D 1 (k q ,k r ,k s )B kq B kr A ks +T\le^ +k ^A kq A kr A ka -T\\e^ +k ^ A kq B k B ks ) = 0, (64) 

q,r,s=l 

£ (D 2 (k q , k r , k s )A kq B kr B ks + Tl \j k °A k A kr A ka - Y\ \e ik ° A kq B kr B ks ) = 0, (65) 

q,r,s=l 

where 

Di(h, k 2 , h) = rl ° + e ^+^a) r o g _ e * ( r o o + r 2 \ - r°\ - r 2 g) , 

£>2(*i, fc, fes) = rg ° + e^ i+fc2+fe )r? J - e^ 2+fc3 ) (r° ° + r} \ - r\ ° - r° 2 ) . (66) 

Differently from the n = 2 case the new relations (65) and (66) are in general not consistent 
with relations (57) and (58). This will restrict the integrability of the Hamiltonian (39) on 
special manifolds of the coupling constants T^ n . 

As in the n = 2 case, the use of (60), with n — 3, imply that an arbitrary amplitude 
of |^3,p) is proportional to Tr A kl A kj A km E L VL P . The cyclic property of the trace supple- 
mented by the algebraic relations (43), (49) and (59) give us the relation that fix, the up to 
now free, complex spectral parameters {k±, k 2 , k 3 }, i. e., 
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e ^ = _JJ 5 (fc i ,fc I ), (j = l,...,n), 



(67) 



with n = 3 and S^/cj, fcj) given by (59). 

n = 4. The relations coming from the amplitudes with collisions with up to three 
particles are solved by (48) and (53), with energy and momentum given by (50). The 
algebraic relations obtained are generalizations of (57), (58), (64) and (65). A new relation 
relating the product of four matrices comes from the amplitude f(x,x,x + l,x + 1) where 
four particles are at colliding positions: 



e 4 Tr 



E X ~ 1 BE~ 1 BBE~ 1 BE L ~ X ~ 1 Q P 



— F 1 x Tr 
i 02 il 



-r 2 °Tr 

i 02 il 



E x - 2 AABE- x BE L - x ~ x ^l P 



— F 1 x Tr 

1 2 11 



E x - x BE- x BAAE L - x - 2 Sl P 



E x - 2 BE- l BEBE~ l BE L - x - x Sl, 



- T° 2 Tr 

1 2 11 



E x ~ 1 BE~ 1 BEBE~ 1 BE L ~ x ~ 2 Qj 



I ( q-pO _ r 2 _ r 2 _ r 2 2\ T 
[° L 00 i 02 i 20 1 2 2,1 11 



E x - x BE~ x BBE- x BE L - x - x Sli 



(68) 



Inserting (48), (50) and (53) and using (57) and (58) in this last expression we obtain the 
algebraic relation with four spectral parameter matrices 

4 



(D 4 (k q ,k r ,k s ,k t ) B kq B kr B ks B kt 

q,r,s,t=l 

+Y\le^ +k ^A kq A kr B ks B kt + ?l\e ikt B kq B kr A ks A kt ) = 



(69) 



where 



D 4 (h, k 2 , fc s , k 4 ) = rg jj + T° g e *<*i+fa+*»+**> _ ^2 2 + r _ r 2 _ T 2 0) ^(fcs+M (7Q) 

This relation will impose, in addition to (64) and (65), further restrictions for the coupling 
constants of the general Hamiltonian (39). The spectral parameters (hi, . . . , fc 4 ), as in 
the previous cases, are fixed by the cyclic property of the trace, and are given by (67) with 
n — 4. 

General n > 4. All the amplitudes produce relations that are solved by the spectral 
parameter matrices introduced in (48) and (53). These matrices {^U.,} and {B kj } (j = 
1, . . . , n) should obey the set of algebraic relations 
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(A kj ) =(B kj ) =0 (j = l,...,n) 



(71) 



Y, N (Ki,KMk P1 A kp2 =Y, C i(^K*) B k P1 B ^ (jy^l = l,...,n) (72) 

p(2) p 2 

]T C (k Pl ,k P2 )B kpi B kp2 = ]T C 2 (k Pl ,k P2 )e ik ^A kpi A kp2 (j ^ 1 = 1, . . . , n) (73) 

n(2) P2 



E 

p(3) 



X! (D^k^kp,, k P3 )B kpi B kp2 A kps + Y\le l{k ^ +k ^ A kn A kv2 A kpz 

p(3) 

_ r i Je^+^J^s^S^) = 0, (74) 

(D 2 (A; P1 , A; P2 , A; P3 ) A kpi B kp2 B kp3 + rj \e*** A kp± A kp2 A kpa - F\ \e ik ™ A kp± B kp2 B kpa ) = 0, 

(75) 

£ (^(Aw, fcp 3 , K,)B kpi B kp2 B kp B kp4 + I* Je'^+^+^U^^fl^B^ 

p( 4 ) 

+r^e^ 5fepijBfcp2 A fep3 A fcp4 )=0, (76) 

where in the above expressions the sums are over the permutations p^ (m = 2,3,4) 
{pi, . . . ,p m } of m distinct integers . . . ,j m ) taken from the set (1, 2, ... , n). Since The 
general Hamiltonian (39) has only nearest-neighbor interactions we have no new relations. 
In fact only (74)-(76) produce constraints for the integrability of (39). Exploring the fact 
that Ci(k, k') is symmetric through the interchange k <-> k' we can use (72) to eliminate the 
matrices {B k } in (74)-(76): 

E {-^i(^pd ^P2! kp;i)N(k Pl , k P2 )Ci(k P2 , k P3 ) + e l<ykp2+kvz ^ T\\ ) Ci{k P2 , k P3 ) — T\\N{k P2 , k P3 ) 

p(3) 

xC^k^kp^dik^k^A^A^A^ =0 (77) 
E {^(^pd k P2 , k P3 )N(kp 2 , kp 3 )Ci(k Pl , k P2 ) + e lfcp3 rJjCi^p!, A; P2 ) — T\\N(k pi , k P2 ) j 

p(3) 

xd(A; P2 , kpjdikp,, k P3 )A kpi A kp2 A kps = (78) 
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{A(^d k P2 , k P3 , k Pi )N(k Pl , k P2 )N(k P3 , k PA ) + r 20 e ^ P2+ P3+ P4 ^Ci(k Pl , k P2 )N{k P3 , k PA ) 

+rj 2e lfep4 N(k Pl , kp 2 )C\{kp 3J fc P4 )| Ci(k Pi: k P3 )Ci(k Pi: k P4 )Ci(k P2 , k P3 )Ci(k P2 , k P4 ) 
xA kpi A kp2 A kp[i A kp4 =0. (79) 

The above relations should be consistent with the commutation relations 

A fc .A fei = S(/^)A fei A fcj (j^Z), { A k 3 ) 2 = (j,Z = l,2,...,ra). (80) 

for any values of kj G C (j = 1, . . . , 4). The use of (80) in (77)-(79) enable us to express the 
left side as a polynomial on the variables e ifcp ? (j = 1, . . . ,4). A sufficient condition for the 
integrability is obtained by requiring that all coefficients of this polynomial are zero. 

Although we did not consider the problem of finding all the possible solutions of (77)- 
(79) with (80) we verified that all the known exact integrable spin-1 chain with a single 
conservation law are solutions of these equations, defining properly an associative algebra. 
These are the cases of the following models. The Fateev-Zamolodchikov model [30] where 

pO _ p2 2 _ n pO 1 _ p2 1 _ p2 _ -nO 2 _ _i pi _ pO 1 _ p2 1 _ pi 2 _ i 
i 00 ~ 1 2 2 — u 5 i 01 _i 21- i 02 _i 20 _ l i 1 1 ~ 1 1 ~ 1 1 2 _ 1 2 1 _ i i 

rhl = r?? = r?? = r|S = 2cos( 7 ), rgl = r 2 ° = -3 + 4sin 2 ( 7 ), 

r^ = r^ = _l + 4sin 2 ( 7 ), r^ = -2 + 4sin 2 ( 7 ), (81) 
and 7 is a free parameter; the Izergin-Korepin model [31] where 

pi _ p0 1 _ pi 2 _ r 2 1 _ 1 p0 _ p0 1 _ pi _ n r 2 _ p0 2 r< sh : " 1 
x 01 _ i 10 ~ i 21 — 1 1 2 — L i ^0 ~ i 01 ~ i 10 ~ U > i 02 _ 1 2 



cosh(37) ' 

t^2 o _ -ni i _ c °s h ( 2 7) ,2 7 t^ii_-pQ2_ cosh(2 7 ) ^_ 27 ±1 _ O cosh( 7 ) cosh(2 7 ) 

J-ll - i 20 — i/Q \ e ' i 02 _ i ll _ i/Q \ e J L ll — Z i/o \ ' 

cosh(37j cosh(3 7 ) cosh(3 7 ) 

p0 2 _ p2 _ p2 , opl 1 ■ I1 i.2 ( ' \ pi 2 _ pO 2 |_ p2 .,-47 p2 1 _ p0 2 , p2 0^47 
1 2 — 1 2 — i 02" t " zi ll blml W )■> 1 12 ~ 1 2 T 1 o 2 e ) i 21 — i 02" t " i 2 e i 

r 22 = 2(r° 2 + cosh(47)r 2 °), (82) 
and 7 is also a free parameter; and the spin-1 Hamiltonian introduced in [32] where 

pO _ p0 1 _ pi _ n pi _ p0 1 _ _i p2 _ p0 2 _ . p2 1 _ pi 2 _ _ 
l00 ~ x 01 ~ 1 1 — u ' 1 01 ~ 1 10 ~ L i 1 02 ~ 1 2 — L i 1 12 ~ 1 21 ~ fc ) 

2 - e 
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rg| = i18 = ~(^ £ + ^), (83) 

with e = ±1 and £ a free parameter. In all these cases the eigenenergies and momentum 
are givem by (50) with {kj} given by (67) with the appropriate S(kj,ki) matrix given in 
(59), in agreement with the Bethe-ansatz solutions of these models. In order to conclude 
this section we mention that analogously to the generalyzed XXZ presented in section 2, we 
can also produce straightforwardly a generalization of the models (81)-(83) where the spin 
1 and spin particles (s z -basis) have now an effective size s (1,2,. .. ) distinctly from the 
size 1 of the particles with spin —1. A choice of the matrices as in (27), i. e., 

n n 

A (s) = J24 S -E 2 - S , £ (s) = Y,B^E 2 -\ (84) 

3=1 3 3=1 

will give the same relations derived previously except that the spectral parameter, instead 
of been fixed by (67) is now given by 

n / ik \ s ~ 1 

^ L = II Stk 5 ,h\ (J = l,...,n). (85) 

1 = 1, Ift \ e / 

IV. MODELS OF SPIN 1 WITH TWO CONSERVATION LAWS. 

In this section we are going to formulate our matrix product ansatz for models describing 
the dynamics of two types of particles on the lattice, where the total number of particles of 
each type is conserved separately. Integrable models on this category are the spin-1 quantum 
chains like the anisotropic Sutherland model [33] or Perk-Schultz model [34] and the t-J 
model [35], and the stochastic Hamiltonian that merge from the problem of asymmetric 
diffusion of particles hierarchically ordered [43] . Similarly as we did in section 1 in order to 
illustrate the versatility of our matrix product ansatz we are going to derive the extensions 
of the above models to the case where the two types of particles (type 1 and 2), or s z = 1 
and s z = in the spin representation of the models have arbitrary hard-core interactions. 
Particles of species 1 (2) will have an effective size s± ($2) that exclude the presence of 
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particles on its rightmost si — 1 sites (s 2 — 1 sites), where si,s 2 = 1,2, The above 

mentioned integrable models correspond to the particular case where Si = s 2 = 1. 

Let us attach to each site i — 1, 2, . . . , L of the lattice, a site variable Qi. that takes 
the values Qi — if the site is empty or the values Qi = 1, Qi — 2 if the site is occupied 
by a particle of type 1 (size si) and 2 (size s 2 ) respectively. The general Hamiltonian we 
consider has an U(1)<E>U(1) symmetry and governs the fluctuations of the configurations 
{Qi, ■ ■ ■ iQl) on a ring of perimeter L: 

KT u(1) = -£^ + Lr™ 



E 3 =V 



2 



.o=l 



E( r oX a <i + r °o<^i) 

V, (86) 



a=l /3=1 a=0 /3=0 



where £ Q,/3 (a,/3 = 0,1,2) are the usual 3x3 Weyl matrix with i,j elements (E ' m ) = 
Si,iS m j and 1^™ are the coupling constants. The projector V in (86) projects out from 
the space of configurations those where the particles occupy forbidden positions due to 
their sizes. The last sum in (86) accounts for the "static" interactions while the first and 
second sums are the "kinetic" terms representing the motion and interchange of particles, 
respectively. The U(1)®U(1) symmetry supplemented by the periodic boundary condition 
of (86) imply that the total number of particles ni, n 2 = 0,1,2,... on class 1 and 2 as well 
the momentum p = M (/ = 0, 1, . . .,L-1) are good quantum numbers. 

We want to formulate a matrix product ansatz for the eigenvectors ^ ni ,n 2 ,p of the eigen- 
value equation 

(87) 

belonging to the eigensector labeled by (ni,n 2 ,P). These eigenvectors are given by 

where the kets \xi, Q±; . . . ; s n , Q n ) denote the configurations with particles of type Qi {Qi = 
1, 2) located at the positions Xi (xi — 1, . . . , L). The total number of particles is n = rii + n 2 . 
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The summation {Q} = {Qi, . . . , Q n } extends over all the permutations of n numbers {1, 2} 
in which n\ terms have value 1 and n 2 terms the value 2, while the summation {x} = 
{x±, . . . ,x n } extends, for each permutation {Q}, into the set of the nondecreasing integers 
satisfying 

> Xi + s Qi = 1, . . . ,n - 1 
sq 1 < x n - x x < L - s Qn . (89) 

The matrix product ansatz we propose asserts that the amplitudes of an arbitrary eigen- 
function (88) are given in terms of traces of the matrix product 

f(x 1: Qi, . . .;x n ,Q n ) = Tr 

The matrices are associated to the particles of type Q (Q=l,2). As in the previous 
sections the matrix E is associated to the vacant sites and the matrix Qp satisfying 

Ett P = e- iP Q P E, Y^ Q) Qp = e - iP Q P Y {Q) (Q = 1, 2) (91) 

ensures (see (7)) the momentum P = -jr-l (I = 0, 1, . . . , L — 1) of the eigenvector. Let us 
consider initially the simpler cases where n = 1 and n = 2. 

n = 1. We have distinct equations depending on the type Q — 1, 2 of the particle. The 
eigenvalue equation (87) give us 

E x Y (Q) E L-x-i np 



_l_ M r 0Q r Q 0\ T 
+ \ zi 00 _i 0Q _i Q0j lr 

where = e^o, = e ,i are the eigenvalues. As in the previous sections a convenient 
solution is obtained by introducing the spectral parameter dependent matrices 

Y {Q) = Y {Q) E 2- SQ (q = i >2 ), (93) 

with k G C, that satisfy the commutation relation with the matrix E 

EY (Q) = e ik Y (Q) E (q = ! 2). (94) 
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^x 1 -ly(Qi) E X2-xi-l Y (Q2) . . . ^n-iTi-l-ly(Qn) E L ~ Xn Vti 



(90) 



£ (Q) Tr EX -1 Y (Q) E L-x n 



-pQ Orp 

-i 0Q lr 



^x-2y (g) E L-x+l fi 



-pO Qrp 

1 Q0 ir 



^-ly (Q) e L-xq 



(92) 



Inserting (93) and (94) into (92) we obtain 



e «> (k) = - (r? Q V* fe + v^ k - 2rg g + rj g + rg J) (q = 1,2) 



27T 

p = k = —l (I = 0,1,. ..,L- 1). 

i7 



(95) 



n = 2. For two particles of types Qi and Q2 (<3i, Q2 = 1, 2) on the lattice we have two 
types of relations. The eigenvalue equation applied to the components where the particles 
of class Qi and Q 2 are at positions (x±, x 2 ) with x 2 > x 1 + sq 1 give us the generalization of 
(92) 

£ (Ql,Q2)rjT r ^Xi-ly(Qi) £iX 2 -Xl-ly(Q 2 ) E L ~ X2 fi,p 
£ixi-2y(Qi) ^X2-xiy(Q 2 ) E L ~ X2 Vlp 
^xi-ly (Qi)^X2-xi-2y (Q2)^i-X2+1^^ 
^1X1— ly(Ql)^rX2— xiy (Q2) jT}L-x 2 -lQp 



-pQl Orp 

_i OQi ir 

pQ 2 Orp 

_i 0Q2 ir 



-pO Qirp 

" 1 Oi o lr 



^xiy(Qi)^X2-xi-2y(Q 2 ) E L ~ X2 Vti 



pU V2rp 

i Q 2 ir 



+ (4r° 



_ r Qi 



-.Qi F 0Q 2 



_ pu i^i pyi u pu 

L 0Qi L QiO i 0^2 



0~j rp^ ^ X1 _ly (Ql)^X2— Xl-ly (Q2) £jL-X2Q^ 



(96) 



while the components where the particles are at the colliding positions (xi = x, x 2 = xi+sq^ 
give us 

£ (Ql,Q2)rp r ^ix-ly(Qi) ^SQ 1 -ly(Q 2 ) £jL-x-s Ql Qj 

-rJ^°Tr r^-syWi) y«a) E L ~ x - s ^n P 



-pO Q2rp 

_i Q2 ir 



£ix-ly(Qi) £S Ql y(Q 2 ) £i£-x-s Qi -1q^ 
£ix-ly(Q 2 ) ^SQ 2 -ly(Qi) E L ~ x - s Q2Q lp 

+ (3r° 8 — rg §i — rg= g — rg; g) Tr [e x ^y ^e s ^ -± y ^e l - x - s ^ q p 



1 Qi Q2 11 



(97) 



Let us consider initially the case where the particles are of same type. In this case we 
have exactly the same situation as in the XXZ chain considered in section 1, and a solution of 
(96)-(97) is obtained by identifying as composed by two spectral parameter dependent 



matrices and Y^\ i. e., 



2 

y«) = ^y«) jE ;2-.g with EY£ Q) = e^Y^ E, (y£ Q) ) =0 (Q = l,2), (98) 
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with n = 2. These last relations, when inserted in (96) give us the energy and momentum 
in terms of the spectral parameters kj (j = 1,2) 



n 



£ (Q,Q) = Y /£ (Q\k J ), P = Y. k i, (99) 

3=1 3=1 

where n = 2 and e^ Q \k) is given by (95). Using (98) and (99) in (97) we obtain the relations 
Y^Yt? ] = S^MYt^ (j^l), (nf) 2 = 0, (l<j,l<n) (100) 

where 

r Q . r Q i(k,+ki) ( pO _i_ p<9 Q p<2 r Q\ ik, 

cQQ (h m_ L QQ + L QQ e ^oo + ^Q + QQ L QQ) e , _ , , . 

*QQ{Kj,Kl)- r Q0 r 0Q i{kj+kl) _ ( r00 + r QQ r Q _ r Q\ ikl W UUlj 

Let us consider now the case where the particles are of distinct species. The distinguibility 
of the particles allows two type of solutions of (96)- (97). We may tray a standard solution 
as in (98) (case a) where each of the matrices (Q — 1,2) are composed by two spectral 
parameter matrices, with the same value of the spectral parameters ki,k 2 (see (98)) or 
alternatively (case b) we may consider a special solution where each is composed by a 
single spectral parameter matrix, with a distinct spectral parameter, i. e., 

= y(V 1 , F (2) = Y k fE 2 - 82 with EY k f = e ik iY k fE (j = l,2). (102) 

In case b (96) give us the energy and momentum as in (99) while (97) give us two independent 
equations: 



P<32 . r 0Qi J(ki+k 2 )\ , (-nOO i pOi Q-2 pQi o r 0Q 2 \A 
i 0Q 2 " hi Qi0 e ) ^ 00 1 Qi Q 2 J^QiO i 0O 2/ l e 



(Qi)vWa) 



2 



+r£ S^^S^ = o (Qi ^ Q 2 = 1,2). (103) 

Since at this level we want to keep k\ and ki as free complex parameters (103) imply special 
choices of the coupling constants Tf™ of (86), for example 3 , 

p2 _ p0 1 _ p2 1 _ n p0 , pi 2 _ pi , p0 2 (-[(\A\ 
i 02 — i 10 ~ i 12 — u ) i 00" t " i 12 — i 10" t " i 2> l iLK V 



3 Another equivalent choices are obtained by the interchange 1 <-> 2 in (104)-(106). 
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that gives 



where 

Sl 2 2 (k u k 2 ) 



v (l) v (2) _ 012/L K \ V (2) V (1) 
v (2) v (l) _ rai/, , n v (1) v (2) 

pi 1 pO 2 j(fci+fe 2 ) _ (pO , p2 1 _ p2 _ pO 1\ ifci 
i 01" l_i 2 e l i 00" l_i 21 i 20 1/ 



(105) 



(106) 



Let us consider case a. A solution of (96) and (97) where (Q = 1, 2) is given in terms of 
two spectral parameter matrices {Y^} as in (98) is possible only if these matrices satisfy 



y(Ql)y(Q 2 ) 



j = 1,2 (g 1 ,g 2 = i,2), 



(107) 



and the Hamiltonian (86) has its coupling constants restricted to 



r l _ r 2 r 1 _ r 2 
i 01 ~ 1 2) 1 10 ~ i 20' 



With (107) and (108) the energy and momentum are given by 

£ (Qi,Q 2 ) = e Wi)( fcl ) + e«»)(jfe 2 ), P = h x + fc 2 , 



(108) 



(109) 



2 

E 

l,m=l 



J- r>. n„C 



1 pikn 
1 V2 



y(Qi)y(Q 2 ) 
y(Q 2 )y(Ql) 



0. 



(110) 



respectively. Inserting (98) and (109) into (97) we obtain algebraic relations that can be 
written in a matrix form 

■nQi Q2 ik m 7-) 1 ifc„ 

where 

A,™ = - (r %° + r° Q V (fei+M ) , v Q2 , Ql = rgg + rg; g - rg; g - r° g (q 1; q 2 = 1,2). 

(in) 

This last relation can be rearranged straightforwardly (see e.g. [43]) giving us (Qi 7^ Q2) 



ta(Qi)^(Q2) _ qQlQ2(u u \y(Q2)y(Ql) 1 cQi Q2/K K \y(Ql)yC?2 
r fe m ~ D Qi Q 2 K mJ r fc m + D Q 2 Qi K rn)Ik m 1 k t 



(112) 



where 
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SS§(*|. U = - I 1 [(^,m + ^ llQ2 e^) ^ 2 , Ql - ge*» 

1 [(2>i,m + v Qi,Q2 eikm ^ * - g^ gilQ2 e^] | (113) 

and 

A = (TV + ^ 2lQl e^) (l> I>m + v QuQ2 e ik -) - V Q Q \ grg ge^. (114) 

In all cases, for Qi = Q2 (see (100) and (101)) or Q x 7^ Q2 in case b (see (104)-(106)) or case 
a (see (112)-(114)) the cyclic property of the traces in our matrix product ansatz (90) will 
fix the spectral parameters. Instead of producing these equations for n = 2 let us consider 
the case of general values of n. 

General n. We now consider the case of arbitrary numbers ni, n<i of particles of type 
1 and 2 [n = n\ + n 2 ). The eigenvalue equation (87) when applied to the components of 
the eigenfunction \^ ni ,n 2 ,p) where all the particles are not at colliding positions, give us 
a generalization of (92) and (96) that can be solved in two distinct ways as in the case 
n — 2. In case b, where the coupling constants satisfy (104), we identify the matrices 
(Q = 1,2) as composed by the spectral dependent matrices 

ni n 

y(i) = ^7«^- Sl) y(2) = Y^E 2 ~ S \ (115) 

j=l j= ni +l 



where Y£ satisfy 

EY (Q) = e *iyW) s> (F^) 2 = 0, hq-! < j < n Q + (Q - l)ng_!, (Q = l,2), (116) 

with n = 0. On the other hand in case a, where the coupling constants satisfy (108), the 
matrices Y^ are given by 

n 

Y (Q) = Y J Y^ ) E 2 ~ S Q, with EYj® = e ikj Y^E 1 < j < n. (117) 
3=1 

In cases a and b the energy and momentum, in terms of the spectral parameters {kj}, are 
given by 
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w 2 = E £Ql (%) + E £ ° 2 (%), ^ = E*i- ( 118 ) 

jr' = l j=ni+l J = l 

The amplitudes of \^ ni ,n 2 ,p) where a pair of particles of types Q\ and Q 2 are located at 
the closest positions Xi and = xt + sq 1 will give us the following algebraic relations. In 
case b we obtain 

Y^Y^ = S^ikj, h)Y™Y™ (Q 1 ,Q 2 = 1,2) (kj ^ k), 

n Ql -i <j< n Ql + (Qi - n 02 _i < Z < n Q2 + (Q 2 - l)ng 2 _i, n = 0, (119) 

where the algebraic structure constants are the diagonal S'-matrix with non-zero elements 
given by 

r Q . r 0,Q A(ki+k,) (-pOO ,-pQQ r Q0 r 0Q\ iL 
qQQ(u u\- 1 QQ + 1 Q,o e ~ { L 00 + L QQ~ L Q0~ L 0Qj e J 

,12. j \ _ 1 _ r o i + r ° ge^ +fc ») - (rgg + rfj - r| g - rg j) 

61 2 - su(ki,kj) - • { } 

In case a, where the coupling constants satisfy (108) we obtain 

^'if^EE^^gif 1 ^ (^U, l</,m<n, (121) 
o;=i q' 2 =i 

where the non vanishing values of the non-diagonal S'-matrix are the six values given in 
(113) and (101) with the choice (104). 

For arbitrary amplitudes we have in our matrix product ansatz (90) a product of n ma- 
trices {Y k ^}. Our ansatz will be valid only if the relations (119) in case b or (121) in case a 
provide a unique relation among these products, otherwise \^ ni ,n 2 ,p) m n °t properly defined. 
This means, for example, that the products • • • Y^Y^Y^ ■ ■ • and • • • Y^Y^Y^ ■ ■ ■ 
should be uniquely related. Since we can relate then either by performing the commu- 
tations in the order 0^7 — > fia^ — > f3^a — > ^(3a or af3^ — > a^f3 — > 7a/? — > the 
structure constants S"y of the algebraic relations (119) and (121) should satisfy 

7,7',7"=1 

E S^fik^k^^ik^^Sj^h,^), (122) 

7,7', 7 "=1 
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for a,a',a",P,P',P" = 1,2. This last constraint is just the Yang-Baxter relation [44,2] of 
the S'-matrix defined in (119) and (121). Actually the condition (122) is enough to ensure 
that any matrix product of spectral matrices {Y^} is uniquely related and it implies the 
associativity of the algebra of the operators {Y^}. 

In case b, the S-matrix given by (120) is diagonal and the Yang-Baxter relation (122) is 
satisfied trivially, so that the only restriction in the coupling constants of the Hamiltonian 
(86) is given by (104). In case a the S'-matrix given in (101) and (113) is non diagonal and 
the Yang-Baxter relation (122) produces strong constraints in the allowed couplings of the 
Hamiltonian (86). 

We did not consider in this paper the problem of finding all the possible coupling con- 
stants r™ ; n in (86) that renders the algebra (121) associative, or equivalently, to find all 
the solutions of the Yang-Baxter relations (122). It is important to stress that our matrix 
product ansatz (90) produce the same relations (121) independently of the hard-core sizes 
Si and s 2 of the particles of type 1 and 2, respectively. Solutions of (122) were presented in 
the literature along the years. The solution where 

rj£ = l r^ = 6 Q cosh( 7 ), 

T a J = sign^ - P) sinh( 7 ) (a^0), (123) 

where 7 is a free parameter and e±, e 2 , e 3 = ±1, corresponds to the anisotropic Perk-Schultz 
model [34]. The isotropic model obtained by setting 7 = and ei,e 2 ,e 3 = ±1 is the SU(3) 
Sutherland model [33]. The solution (123) with e\ = — e 2 = €3 = 1 give us the anisotropic 
supersymmetric t-J model 4 . A solution of (122) is also known [47] for the case where 

r?S = -r°? = g si ^ a -v (a ^(3), r^ = o («,/? = 0,1, 2), (124) 

with q real and q > 1. The Hamiltonian (86) with the couplings (124) describes the time 



4 A Jordan- Wigner fermionization of the Hamiltonian (86) with s\ 
(123) give us the anisotropic supersymmetric t-J model [35]. 
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= s 2 = 1 and coupling constants 



fluctuations in the asymmetric diffusion problem of two species of particles hierarchically or- 
dered (asymmetry q). A generalization of this problem, treated through the matrix product 
ansatz is presented elsewhere [45]. 

In order to complete our solutions through the matrix product ansatz (90) we should 
fix the spectral parameters (ki, . . . , k n ). Using the algebraic relations (116) or (117) an 
arbitrary amplitude is proportional to 



Tr 



ki 

The cyclic property of the trace fix the spectral parameters. 



(125) 



Let us consider case b, where the coupling constants satisfy (104). The commutation 
relations (91), (116) and (119) give us 



/ ikj\ si" 1 ' 



n sll^M 

l=ni+l \ c / 



n siK^^r— 



e ik 3 \ «2-l" 



l=m+i 



ni 



n^i(^*i)hib 



e ik 3 \ 



1 < j < ni) 



(m<j<n), (126) 



i=i 



where the S'-matrix is given by (120). The energy and momentum are obtained by inserting 
the solutions of (126) into (118). The Hamiltonian (86) that is exactly integrable in case b 
(see (104)) is given by 

L 

n - 2^' \ L 1^j - C/ i+l + i 20- C/ j ^j+l + 1 2 l^j j+S2 j+si ' ^ 1 l^j ^j+S! 
3=1 

| -p2 2 7Tl2 2 7Tl2 2 | "pO 1 TTlO TTll 1 I pi TTll 1 7710 , pO 2 T7>0 77>2 2 , p2 7T>2 2 TTlO 

+ 1 2 2^j L j+s 2 + 1 A' + 1 1 O^j h j+ S1 + 1 2^j &j+l + 1 2 ^j-M 2 

+r1 1^ 2 e]1 2 + (rj ; + r° g) ^ ^ + 2 S1 } p, (127) 

where we chose Tq [j = but still we have 10 free parameters! The particular choice Si = 

„ _ 1 F 1 2 — pi 2 _ pi , pO 2 pi _ _pl p0 2 _ _p0 2 pi 1 _ p2 2 _ pO 1 _ p2 _ n 

52 — 1, 1 2 i — i 12 — i 01"r i 20> i 01 — 1 1 0' 1 2 — 1 02) 1 1 1 _ 1 2 2 _ 1 1 _ 1 2 _ u i 

gives the Hamiltonian related to the stochastic problem of fully asymmetric diffusion of two 
kinds of particles, whose exact integrability was obtained in [26] through the dynamical 
matrix product ansatz. 

Let us return to the general case. The cyclic property of the trace in (125) and the 
commutation relations (91), (117) and (121) give us 
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Tr [y^ 1 ** . . . yW'O^-k— i)— "2(32— 



gife ; ;[L-ni(si— l)-n2(s2-l)] 



x £ (Qi, . . . , QnlTlQi, . . . , Q^Tr [y fc f^ • • • y«'-) S W(.i-i)-n a («-i) np ] > (128 ) 
Qi,...,0;=i 



where we have used the identity (see (101) and (113)) 



E && 1 (k j ,k s ) = -l, (129) 



and 

(Q 1 ,...,Q n \T\Q' 1 ,...,Q' n ) 

= E {^q; gf ' ' ' Q'f 1 ' ' ' (^«' kj) elP ^ Sl ^ } (130) 

is the transfer matrix of a 2™ x 2 n -dimensional transfer matrix of an inhomogeneous 6- vertex 
model (inhomogeneities {h}) with Boltzmann weights given by (101) and (113). The 6- 
vertex model is defined on a cylinder of transversal perimeter n with a seam along its axis 
producing the twisted boundary conditions 

S Q Q 7 Q Q t{kn,k) = S%l Q Q l(k n} k)<P(s QI( ), (131) 

where 

^(s) = e iP{s - 1 \ (132) 

and, as before, P is the momentum of the eigenstate (88). The relation (128) give us the 
constraints for the spectral parameters 



e 



-^(L+n-n 1S1 -n 2S2 ) = A ^ (j = 1, . . . (133) 



where A(kj, {h}) are the eigenvalues of the transfer matrix (130). The condition (133) leads 
to the problem of evaluation the eigenvalues of the inhomogeneous transfer matrix (130). 
This can be done through the algebraic Bethe ansatz [46] or the coordinate Bethe ansatz 
(see [43] and [47] for example), and we obtain 

A(%, {ki}) = 0( S2 ) n SH{k h kj) n Q \Z 3 , (134) 
i=i i=i^2i(kj,k l ) 
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where the unknown parameters k^ (I — 1, . . . , n\) are fixed by the n\ coupled equations 
<H g i) i~r Sjijkh kj ^) _ -j-j- S%%(kj \ k\ ^) S%l(kl \ kj ^) 

for j = 1, . . . ,tii. The equations (133)-(135) fix the spectral parameters {kj]j = 1, ... ,71} 
of the matrices introduced in our matrix product ansatz. 



V. MODELS OF SPIN-| WITH TWO CONSERVATION LAWS. 

Like in the previous section we now formulate our matrix product ansatz for U(1)®U(1) 
symmetric models describing the dynamics of two types of particles on the lattice, where 
the total number of particles of each type is conserved separately. However now, distinctly 
from the last section, particles of distinct types may occupy the same site on the lattice 
(double occupancy). Models on this category are the spin-| anisotropic Perck-Schultz model 
[34], the Essler-Korepin-Schoutens model [36], the Hubbard model [37], the Hamiltonian 
derived in [48] from the i?-matrix introduced in [49], and the two-parameter integrable 
model introduced in [38]. 

In order to define the Hamiltonians for these models on a lattice with L sites let us attach 
to each lattice site a variable Qi (i = 1, . . . , L), that takes the values Qi — if the site is 
empty, Qi = 1, 2 if the site is occupied by a particle of type 1 or 2, respectively, and Qi = 3 
if the site have a double occupancy of particles of type 1 and 2. In the fermionic version of 
the models the particles of types 1 and 2 are the electrons with distinct spin polarizations. 
The most general model with nearest neighbor interactions and periodic boundary condition 
is given by 

a^/3=0 77^=0 a=0/3=l 

where the coupling constants satisfy T"^ — Oiia + P^j + v, r^" = r^J = unless 
a — P — 7, and E a/3 (a, j3 = 0,1,2,3) are the 4x4 Weyl matrices with i,j elements 
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(E a/3 ) . = SajSpj. Without loosing any generality we also chose hereafter Tqq = 0. The 
first and second summation in (136) account for the kinetic and static interactions. The 
U(l)<g)U(l) symmetry and the translation symmetry, due to the periodic boundary condition 
of (136), imply that the total number of particles ni, n 2 (0, 1, . . .) of particles of types 1 and 
2 and the momentum p = M (/ = 0, 1, . . . , L — 1) are conserved separately. 

We formulate a matrix product ansatz for the eigenvectors \^ ni ,n 2 ,p) of the eigenvalue 
equation 

H uwm ^ m , n2 ,p) = e niina \* nitna , P ) (137) 
in the eigensector (ni,n 2 , P). An arbitrary eigenfunction is given by 

Qn) Ql j • • • j x m Qn) (138) 

{Q} {*} 

where \x±, Q±; . . . ; x n , Q n ) are the configurations whose particles of type Qi — 1, 2 are located 
at Xi = 1, . . . , L and n = rii + n 2 . The summation {Q} = {Qi, . . . , Q n } extends over all 
permutations of n numbers {1, 2} in which ri\ terms have the value 1 and n 2 terms the value 
2. The summation {x} = {x x , . . . ,x n } extends, for each permutation {Q}, into the set of 
non- decreasing integers satisfying 

(x i+1 -Xi)>l if Qi+i = Qi 

(x i+ i-Xi)>0 if Q i+1 ^ Qi, (i = 1, . . . ,n - 1). (139) 

In order to formulate a matrix product ansatz we associate to the sites occupied by Qi = 
0, 1 and 2 the matrices E, and Y^ 2 \ respectively. In analogy to the results of section 3 
we associate to the double occupied sites (Qi = 3) the matrix product = B^E^ l B^> . 
Our ansatz asserts that the amplitudes corresponding to the configurations with no double 
occupied sites are given by the traces 

f(x u Qi, . . .;x n ,Q n ) = Tr 

while the amplitudes related with configurations with double occupancy at Xj+i = Xi are 
given by the traces 
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(140) 



f(xi, Qi, . . . ; Xi, 1; Xj, 2; . . . ; x n , Q n ) — Tr E 1 y^ 1 ) 

. . . ^- x i-^- 1 E^ 1 E^ 1 '^ 1 ■ ■ ■ E L ~ Xn Vlp . (141) 

The general cases are given by generalizations of (140) and (141). 
Similarly as in the previous sections the matrix Q P satisfying 

EQ P = e- lP Q P E, Y iQ) n P = e- lP Q P Y (Q \ B (Q) Q P = e- iP Q P B (Q) (Q = 1, 2) (142) 

ensures (see (7)) the momentum P of the eigenvector |\l/ niin2) p). 

As in the previous sections in order to satisfy the eigenvalue equation (137) we identify 
and B^ (Q = 1,2) as composed by n spectral parameter dependent matrices 5 



Y (Q) = ^ y<Q) Ei B (Q) = ^ B (Q) E(Q = 1? 2)> 

3=1 3=1 



(143) 



satisfying the commutation relations 



K W) £ = e ifc ' EYu® , B ( k Q) E = e ik i EBP , 



yW) n 



y(<3)yW) 



B (Q) B (Q') 

kj kj 



(Q,Q' = 1,2), 



(144) 



where kj (j = 1, . . . ,n) are unknown complex spectral parameters. 
Let us consider initially the sectors with small values of n. 

n = 1. We have distinct equations depending on the type Q — 1, 2 of the particle. The 
eigenvalue equation (137) give us, for Q — 1,2, 



£ (Q) Tl l E x-l Y (Q) E L-x Q 



E x-2 Y (Q) E L-x+l Qp ^ _ r Q Tr 



^y(Q) E L - x - 1 n P 
(145) 



Inserting (143) with n — 1 and using (144) we obtain the energy and momentum as a 
function of the spectral parameters 



5 Similarly as the solutions (104) in section 4 there exist special solutions, that for brevity we do 
not consider here, where the matrices B^ and Y^ , B^ are composed by the distinct sets 
of spectral parameters {ki, . . . , k ni } and {k ni , . . . , fc„ 1+n2 }, respectively. 
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e w)(jb) = - (r? °e~* k + r° g ? e * + r° g + r§ J) (q = i, 2), 



27T 

P = k = —l (I = 0,1,..., L- 1) 

Lj 



(146) 



n = 2. For two particles of type Q\, and Q2 (Qi,(?2 = 1,2) located at x\ and x<i the 
eigenvalue equation (137) produce distinct relations. The components where the particles 
are not at "colliding" positions, x i+ i > + 1, give us the relation 



£ (Qi,Q2)']T r ^xi-ly(Qi)^X2-xi-ly(Q2) E L ~ X2 Qi 



pQl Orp 

_i OQi ir 



pO Qlrp 

1 Qi ir 



^xiy(Ql)^X2-a;i-2y(Q2) E L ~ X2 Vti 



pQ2 Orp 

_i 0Q2 ir 



pU V2rp 

_i Q2 ir 



£ixi-2y(Qi)£.x 2 -xiy(Q 2 ) E L ~ X2 Vtp 
jjjxi—lY (Qi) — xi— 2y (Q2)^£<-X2+1q^_ i 
^1x1 - 1 y (Qi ) ^710:2 -xi y (Q 2 ) J^L~x 2 - 1 ^ 

" ( r o q\ + r S o + r o S + r S o) Tr [£ Xl_1 y (0l) £ X2 - Xl - V 02 ^-* 2 ^] . (147) 

Inserting (140) with n = 2 into this last expression and using (144) we obtain a solution 
provide the coupling constants in (136) satisfy the constraint 



r° 1 — r° 2 pi o _ r 2 o 

1 1 — 1 2 0' 1 o 1 — 1 o 2- 



(148) 



The energy and momentum in terms of the spectral parameters k\ and ki are given by 



£ (QuQ 2 ) = e Wi)(j bl ) + gW^jfej), P = kl + 



(149) 



where e i ~ Q \k) is given by (146). 

In order to proceed it is better to consider separately the cases where Q\ = Qi and 
Qi 7^ Q2- If <5i = Q2, (137) applied to the amplitudes related to the configurations where 
X\ = x, x-2 = x + 1 give us the relation 



£ {Q,Q)^ X [£ , x-ly(Q)y(Q)^L-x-lQ 



PV Urp 

- _i 0Q lr 



E x-2y(Q) EY (Q) E 1 -*- 1 ^, 



-T° Q Q Tr [E x - 1 Y^EY^E L - x - 2 n P ] - (r°g + T%° + r§ §) Tr [£*-iyW)y«J)^ifi p 

(150) 

Inserting (143) and (149) in (150) and using (144) we obtain the algebraic relations among 
the spectral parameter matrices 
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Y^Y^ = S%%{k 3 MYPYff\ (Q = l,2), 



(151) 



where 



oiki 



r oQ e i( kj+kl ) + r o o + ( r ° g + rg £ - rg g) e 

If Qi 7^ Q2 distinct relations merge from the configurations where X\ = x 2 
Xi — x, x 2 — x + 1. For x\ = x 2 = x we have 



(152) 
x and 



£ (i,2) Tr \e x - 1 B^E- 1 B^E l - x VL p 

- r§ jTr Le^y^yw^-*- 1 ^ 



r£ 2 Tr [s^yMy^ 1 --^, 



F 21 Tr 
1 3 11 



£>x-2y(2)y(l) fiL-xQ 



T 12 Tr 

1 3 11 



£ x-ly(l)y(2) E L-,-l Q 1 - r 3 °Tr f E x ~ 2 E~ l E l ~ x+1 Q 



-r° 3 Tr E x B^E- 1 B^E L ~ x ~ l n P — (rg 1 + r§ g) Tr E x ~ l B^ E~ x B^ E L ~ X VL P , (153) 



while for xi = x, X2 = x + 1 we have for Qi 7^ Q2 (Qi, Q2 = 1, 2) 



£ (Qi,Q2) r ]] r [_£ , x-ly(Qi)y(Q 2 ) E L ~ x ~ l Vli 



_ pQl Orrv. 



1 Q2 ir 



^x-ly(Qi) _£y(<? 2 ) E L ~ x ~ 2 Vt) 



F 03 Tr 
1 Qi Q2 ir 



E x B (l) E -l B (2) E L~x-l n 



E x~2 Y (Qi) E y(Qz) E L - X - 1 Q P 

E x - 1 B { ^E- 1 B^E L ~ x ^lp 
g lr Tr r^-iy(Q2)y(Qi)^^-i^ 



F 30 Tr 
1 Qi Q2 11 



(154) 



Inserting (143) and (149) in (153) and (154) and using (144) we obtain the algebraic relations 

E c (kj,ki)B^ b^ = - £ e^^cl^^oy^V^ + ^^.^y^yW], (155) 



3,1=1 j# 3,1=1 3+1 

+rS&e tt «y« a) y« l) }, g^Q 2 (Qi = i,2), 



y(Ql)y(Q 2 ) 



(156) 



where Co, C^, C 2 , C\, C 2 and L> are the symmetric functions 



C (kj, h) = -Tl\(e ik > + e lkl ) - T%e l ^ +kl) (e 1 ^ + e ikl ) 
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_ (-pi i pO 1 i p2 , -nO 2 _ r O 3 _ p3 0\ „i(kj+ki) _i_ p 3 _l_ F° 3„i2(fc J +fc i ) 
\ i 10" t " i 01" t " i 20" , " i 02 ^3 1 3 ,1 e ^ 1 3 1 3 O e ' 

Ci(kj, ki) = -T\l-T\\e^\ C' 1 {k 3 M) = Tll + Tlle^\ 
C 2 (k v k t ) = -Tl\-Tl\e^\ C' 2 (k v k l ) = TH + Tlle^\ 



(157) 



Multiplying the two relations in (156) by Co(kj, ki) and using (155) we obtain the algebraic 
relations expressed in the matrix notation 



2 

E 



7j,i hi 



y(l)y(2) 
y(2)y(l) 



0. 



where 



-Ci(kj, ki)C[(kj, k)e ik < - C (kj, h) [D(kj, h) - {p\ Q + T° Q \ - T\ fj e ik > 



- C 2 (kj, ki)C[(kj, k{) - C (kj, k{)Y\ \ 

-c 2 (k j: h)c' 2 {k v hy k > - c (%, h) [D{k v hi) - (r° \ + v\ ° - r 2 , *) . (159) 



ife; 



This last equation can be rearranged straightforwardly, giving us 



\r{Ql)\r{Q2) _ 0Q1Q2/7, J. \ V^(Q2 ) y (Ql ) , oQi Q2 I J. h, ^yWDyC 
Y ki r k m - &Qi Q 2 \ K h K m)I km 1 k l -T a Q2Qi\ K h K m) 1 k m 1 h 



where 



^2 fcj) 
S\l{kj, ki) 

S\\{kj, ki) 
Multiplying (156) with Q x 



= - 1 + 



. L + a u (Jkj - 6 hi) - jhi (Pij-PuY 



(158) 



(160) 



(161) 



1 by Ci and (155) by Y\\ and subtracting the obtained 
expressions we obtain a relation that express the spectral parameter matrices in terms 
ofvfU e., 
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where 



W = ^^« CJVO, (162) 



L> 3 (%, fcj) = fc,) - (rig - r} 2) ^(fc., h) - C'^kj, ki)Tl \e ik \ 

P {k 3 M) = -Pj,ie- ikl = c^hjM&ikjM) - c (%,^)r?i (163) 



Since the matrices B^j' only appear in the pairs B\.' B^.' in our matrix product ansatz (140)- 
(141) the relation (162) enable us to express the amplitudes only in terms of the spectral 
parameter matrices Y^p (Q = 1,2). Consequently (160), (161) together with (144) give us 
the complete algebraic relations for the matrices appearing in our matrix product ansatz in 
the sector with n = 2 particles. 

The spectral parameters k\ and fc 2 , as in the previous sections, are fixed by the cyclic 
property of the trace defining the amplitudes. Instead of producing the equations fixing 
{kj} for n = 2 let us consider the cases of general values of n. 

General n. We now consider the case of an arbitrary number n 1: n 2 of particles of type 
1 and 2 (n = n\ + n 2 ). The eigenvalue equation (137) when applied to the components of 
\^m,n 2 ,p) with no particles at "colliding" positions (x i+ i > x^ + 1, i — 1, . . . , n) give us a 
generalization of (147) that is solved by giving us the energy and momentum as a function 
of the spectral parameters 



e nun2 = E^(%) + E c (2) (k 3 ) = E (riSe-* + ^h***) 
j=i j= ni +i j=i 

n 

+^i(r^ + r°l)+n 2 (r^ + r^), p = $>. (164) 

3=1 

The amplitudes of \^ ni ,n 2 ,p) where a pair of particles of type Q 1 and Q 2 are located 
at the colliding positions Xj+i = Xi or = Xi + 1 give us algebraic relations that are 
generalizations of (151) and (160), 

nfX Q2) = £ ^S^Y^Y™ (k^ kl ), 
Y k f l] Y^=0 (j,/ = l,...,n), (165) 
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where the nonzero elements of the S-matrix in (165) are those given in (152) and (161). 
Moreover any pair of matrices B^B^ (j ^ I; j,l = l,...,n) are expressed in terms of 
b y (162) 

For arbitrary amplitudes the matrix product ansatz (140) and (141), after using (162), 
contains a product of n spectral matrices {Y^}. The ansatz will be valid only if the 
algebraic relations among these matrices give us a unique relation among these products. 
This imply that products like • • • Y^ l] Y^ 2) Y^ 3) • • • and • • • Y^Y^Y^ ■ ■ ■ should be 
uniquely related. Similarly as discussed in section 4 these products can be related by apply- 
ing the commutation relations in distinct ways. Performing the commutation relations in the 
order Q1Q2Q3 -> Q2Q1Q3 -> Q2Q3Q1 -> Q3Q2Q1 or Q1Q2Q3 -> Q1Q3Q2 -> Q3Q1Q2 
Q3Q2Q1 will impose constraints among the structure constants Sq^ q 2 of the algebraic rela- 
tions (152) and (161), namely, 

Q,Q',Q"=l 

= E S$^}(k 2 ,k 3 )s2l2l(ki,ks)s2 2 %(k 1 ,k 2 ), (166) 

Q,Q',Q"=1 

for Qi, Q[, Q'l, Q 2 , Q' 2 , Q' 2 f = 1, 2. These constraints are just the Yang-Baxter relations [29,2] 
of the S'-matrix defined in (152) and (161). The Yang-Baxter relations (166) imply the 
associativity of the algebra of the matrices {Y^}. These relations will impose severe 
constraints among the coupling constants T^™ of our general Hamiltonian (136). However 
the Yang-Baxter relations (166) are not enough to ensure that the eigenfunctions \^ni,n 2 ,p) 
are given by our matrix product ansatz, since new relations among the spectral parameters 
matrices {Y^} happens when we have on the lattice 3 or 4 particles at matching conditions. 
These new relations happen because within the range of the nearest neighbor interactions 
of the Hamiltonian (136) we may have up to 4 particles. The eigenvalue equation (137) 
applied to the amplitudes where we have a particle Q at x\ = x and a pair at x 2 — x + 1, 
and to the amplitude with a particle Q at x\ = x + 1 and a pair at x give us, respectively, 
the equations 
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( e w) + e (i) + e (2) + r ° g + rg 3 + r 3 °) Tr ^-lyfOBWr^^^-^p 

T^Tr \e x - 2 Y^EB^E- 1 B^E l - x - 1 Q p ] - rf Q Tr [^-iy(Q)y«')y(Q)^-*- 2 flp 



0Q X1 [-^ J ' ^ ^ w,» 



3 



1 3 11 



^-iy(Q) ^5(1) E -i B (2) E L - x - 2 tt P 

, Q^Q' (Q,Q' = 1,2), (167) 



( e «) +e (l) +e P) +r 03 + r 3g + r O0 )Tl . 



1 3 ir 



03 Ti 3O Ti Q0 

E x ~ l B (1) E- 1 B (2) EY {Q) E L - x - 2 tt P 



E x ~ l B {1) E- l B {2) Y iQ) E L ~ x ^ 1 n f 



-r^Tr E^Y^B^E^B^E^^Qp , Q ^ Q' (Q,Q' = 1,2). 



if 3 Q Tr £«-2y«')yW)y«) ^-^Qp 
E x ~ 2 B w E~ l B^ EY {Q) E^^Vtp 

(168) 



F 3 °Tr 
1 3 11 



Inserting (143) and (164) into these last equations, and using the generalizations of (155) to 
simplify the obtained expressions we obtain the four algebraic relations 

E {D^\k 3 ,k h k m )Y^B^B^ 

j,l,m=l 

_ e ik m ( r 3Q B W B (2) Y (Q) _pQ' Qy(Q)y(Q')y(Q) _ pQ Q'y(Q)y(Q)y(Q') > \ I = q 

E {^ (Q ' 3) (%,^,^)« ) n ( f ) 

j,l,m=l 

_ e ^ + M (r? 3y«>s« s£> - if SyWyWy™ - rf ?y™y™y™)} = o, (169) 

where Q, Q' = 1, 2 or Q, Q' = 21, and 

z^) ( ^, fc„ u = i* <g e <fc+*+*»o + (rg° + r° 3 - rg J) e *(*»+*») + r 3 °, 
dw- 3 ^-, fc„ u = r°g e «*+*+*»o + (r°g + r 3 ° - r*g) e "~ + r? Q °. (170) 

In the case of n = 4 particles the new relation comes from the amplitudes corresponding to 
the configurations where we have two pairs of particles at sites x and x + 1. The eigenvalue 
equation (137), when applied to these amplitudes give us 



(2e« + 2e^ + T° 3 + T 3 3 + P 3 °) Tr [e^B^E^B^ B^E^B^E^^Qp 
-T 2 \Tr \E x - 2 Y^Y^B^E- 1 B^E L - x - 1 n I 
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T 1 2 Tr 

1 3 11 



-F 1 2 TY 
1 3 11 



-F 2 lr TY 
1 3 11 



-F 3 °Tr 
1 3 11 



-F° 3 Tr 
1 3 11 



E x-2 Y (1) Y {2) B 0) E -l B (2) E L-x-l Qp 
E x-1 B (l) E -l B (2) Y (1) Y (2) E L-x-2 Qp 
E x-1 B (l) E -l B (2) Y (2) Y (1) E L-x-2 Qp 
EX ^ B (l) E -l B (2) EB (l) E -l B (2) E L-x-l Qi 
E x-1 B (1) E -1 B (2) EB (l) E -l B {2) E L-x-2 n 



(171) 



Inserting (143) and (164) into this last relation and using the generalization of (155) to 
simplify the resulting expression we obtain 

4 



(i) 
h 



E K' 3) (^, ku k m , k )B^B%B^B { 2 + e^+"°) (i*^ 

j,l,m,o=l 

+niY^>) b«s« + «*Bg) B » (r^y£>y« + r;|y«y( 2 ')} = o, (172) 

where 

D^ 3 \kj, k, k m , k ) = r 3 ° + r° 3 e <(*i+*«+*™+*o) + (r 3 ° + r° 3 - r 3 3 ) e i{ - km+k °\ (173) 

Since any pair B^B^ {k 3 ^ ki) is expressed in terms of the pair Y^Y^\ through (162), 
we rewrite the expressions (169), (170) and (171) only in terms of the matrices {Y^}. 
Multiplying (169) by the symmetric function p(kj, ki)p(ki, k m )p(k m , kj) we obtain 

£ {D^ Q \k 3 , k h k m )D 3 (k h k m ) P {k 3 , h)Ytyt?yt 2 2 

j,l,m=l 

-e^Y%%D,{k 3l kMki, k m )YtfYty^ 

+ (If ?Y™Y™Y™ + r? fYWY^YgP) e ^ p{kl , k{)p{k h k m )} p(k m , k 3 ) = 0, 
£ {D^ 3 \k 3 , k h k m )D 3 (k 3 , h)p{k h k m )Y^Y^Y^ 

j,l,m=l 

-e^'l?^ k m )p(k 3 , h)Y^Y^Yt 2 2 

+ (if ?r fc fV^if) + 1? VYWYgnyW) e^ P (k 3 , kMk, U } p(^, *,) = o, 

(174) 

where (Q,Q') = (1,2) or (2,1). Multiplying (172) by the symmetric combination 
k t )p(ki, k )p(k , k m )p(k m , k 3 )p(k h k m )p(k 3 , k Q ) we obtain 
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£ { [(^ (3 ' 3) (kj, k h k m , k )D 3 (kj, h) + Y\ le^ +k - +k ^p{k v fc,)) L> 3 (^, 

j,l,m,o=l 

+T 2 l e ^m +ko ) Dz{km ko)p{kv kl )YtfYtyty^ 

+Tl\e ik °Dz(k j ,k l )p{k m ,k )Y^ = 0. (175) 

In order to have a matrix product ansatz for the Hamiltonian (136) the relation (174)- 
(175) should be consistent with the two words commutation relations (165), for any values 
of kj E C (j = 1, . . . ,4). The successive use of (165) in (174)-(175) allow us to rewrite the 
left side of these equations as a polynomial on the variables e lkj (j = 1, . . . , 4). Since we do 
not want, on this level, to fix the spectral parameters {kj}, we should impose that all the 
coefficients of these polynomials are zero. This will give further constraints on the coupling 
constant T^™ besides (148) and those imposed by the Yang-Baxter relations (166). 

It will not occur any additional constraint for n > 4. Although we did not consider 
here the problem of searching all the possible solutions of the general Hamiltonian satisfying 
(166) and (174)-(175), we verified that those equations are satisfied by several know exact 
integrable chains. The solution where 

r /j £ = !> r a % = sig n (« - P) sinh( 7 ) - e cosh( 7 ) (a ^ f3) 

rS2 = (ea-eo)cosh( 7 ), (a, f3 = 0, 1, 2, 3), (176) 

where €0,61,62,63 = ±1 and 7 a free complex parameter give us the anisotropic spin-| 
Perk-Schultz model. The fermionic version, obtained by a Jordan- Wigner transformation, 
of the Hamiltonian (136) with the choices (176) and with e = — e 1 = — e 2 = e 3 = 1 is the 
anisotropic version of the Essler-Korepin-Schoutens model. 
The solution where 

r§8 = rg§ = o, rsi = r|S = -¥ (i + M, (177) 
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with U a free parameter, give us after a Jordan- Wigner transformation, the standard Hub- 
bard model with Coulomb on site interaction U [37]. The solutions where the non-zero 
coupling constants are 



r« _ rO a _ i (n^zO) F 3 1 — F 1 3 — F 



3 2 
2 3 



p2 3 
1 3 2 



F3 i 
e ) 1 3 1 



p3 i rl2 



3 



_l_ F 

+ 1 1 2? 



r?^ = = e r^ = e r° 3 = rj§ = r? 3 = er 3 ? = c r§5 = sin(0) 

Tjl = Til = ~er° 2 1 = -erg ° = -|r 3 § = r 2 ^ = rl 2 e" 2 ^ = cos(^) 



p3 2 _ p3 _i_ p2 1 pi 3 _ pi 2 p2 3 _ p2 1 p3 3 _ p3 _i_ p2 1 _i_ pi 
3 2 3 2 1' i 13 ~ i 12> i 23 — i 2l5 1 3 3 — 1 3 1 2 1 "r 1 1 



(178) 



and 



r«0-r0«-i ( n + n) p3 1 _ r i 3 _ F 3 2 _ F 2 3 _ 

L 0a — L aQ — 1 \ a T u )i 1 13 _ ^31 ~ i 23 — i 32 _ e ; 



r? I = T\ I = eTl \ = eTl 3 = ee 2 T 2 = ee 2 T? 3 = e^Tl ? = e^Tj J = sin(0) 
Tl I = rl l = -eT° 2 = - e r 2 § = T 2 = r\ \e~^ = co S (6) 

p3 _ OF2 , s i n 2//,\ (e^-ee-^)_ p 3 1 _ -p3 2 _ -p3 , pi 2 
1 3 ~ Zi 2 + Sill 1^(7; cos( -0) , J-31 - 1 3 2 — 1 3 ^ 1 1 2> 



1 3 
3 



cos(e) 

-i2 1 r3 3 



p2 3 _ p2 1 p3 3 _ p3 U _i_ p2 1 _i_ pi 2 
i 23 — i 2 1) 3 3 3 2 1 ' 12' 



(179) 



where e = ±1 and 9, rj are free parameters, give us the two-parameter integrable models 
introduced in [38] . It is interesting to mention that these two models contain as special cases 
the Hubbard model, the Essler-Korepin-Schoutens model, as well as the g-deformation of 
the extended Hubbard models introduced in [48,50]. 

In all the cases the spectral parameters (fci, . . . , k n ) are going to be fixed by the cyclic 
property of the trace. The relations (155), (151), (162) and (165) imply that any ampli- 



tude of the matrix product ansatz (140)-(141) are proportional to Tr 



Y 



(Qi) 



fci 



■■■Yp n) E L Q 



Successive applications of the commutation relations (151), (165) and (142) give us 



Tr 
xTr 



= e 



ikjL 



]T (Q u ...,Q n \T\Q' v ...,Q' n ) 



Y 



(Qi) 



fci 



Yt Q ' n) E L Q P 



(180) 



where we have used the identity 
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(181) 



and 



(Qu...,Q n \T\Q' 1 ,...,Q' n ) 




(182) 



Q 



is the 2 n x 2™-dimensional transfer matrix of an inhomogeneous 6-vertex model (inhomo- 
geneities {kj} along the vertical) with Boltzmann weights given by (152) and (161). The 
vertex model is defined on cylinder of transversal perimeter n and periodic boundary con- 
dition 



The relation (180) fix the values of the spectral parameters {kj} as the solution of the 
equation 



where A(kj, {h}) are the eigenvalues of the transfer matrix (182). In order to complete our 
solution we need to evaluate the eigenvalue of the inhomogeneous transfer matrix (182). We 
are in a similar situation as in section 4, except that now the inhomogeneous vertex model 
is defined on a periodic lattice. The eigenvalues A(kj, {h}) is obtained by setting in (134) 
and (135) (f>(s) = 1. The spectral parameters k±, . . . ,k n are then fixed by the equation 



Q'nQ'A ^ ' — ^Q' n Q'^ \Kni K)- 



(183) 



MkjAki}) (j = l,...,n), 



(184) 



i=i i=i 



n^2 2 (*i,*i)n 



S^jjkj, k\ •*) 
Sf i(%, k\ -*) 



(185) 




(186) 



for j = l,...,m. 
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VI. THE GENERALIZED XXZ CHAIN IN AN OPEN CHAIN 



Differently from the former sections where the quantum chains are all defined on periodic 
lattices we consider here the formulation of the matrix product ansatz in an open chain. We 
consider the generalized XXZ chain defined in section 2 in an open lattice with diagonal 
z-magnetic fields acting only at the surface points. The Hamiltonian we want to solve is 
given by 

H s = -V s ji £ 0,X +1 + <^<-i + + + h R <r z L }v s , (187) 

where a x , a y , a z are spin-| Pauli matrices, A is the anisotropy and h^, h R are the magnetic 
fields acting at the end points 1 and L, respectively. The projector V s , as in section 2, 
projects out from the associated Hilbert space the configurations where any two up spins 
are at distances smaller than s (s=l,2,. . . ). The choice hi = h L = corresponds to the free 
boundary case. 

The exact solution of (187) for the case s = 1 was obtained through the coordinate 
Bethe ansatz in [39] and through the quantum inverse scattering method in [51]. It is also 
interesting to mention that, in the case where s = 1 the Hamiltonian (187) is SU (2 ) 9 - invariant 
if the anisotropy and surface fields are related by A = (q + l/q)/2, h L = q, h R = l/q. 

Since the Hamiltonian (187) commutes with the total spin operator S z = Ya=\^1 the 
number of up spins n is a good quantum number. We want to solve the eigenvalue equation 

H\V n ) =e„|*„) (188) 

where 

l*n) = • • - . . .,x n ). (189) 

Here x 1 , . . . , x n denote the configurations of the up spins on the chain, and the summation 
extends over all sets of n increasing integers satisfying 

xi > 1, x n < L, x i+ i>Xi + s, (i — 1, . . . ,n — 1). (190) 
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As in section 2 in order to formulate our matrix product ansatz we associate the matrices E 
and A to the sites occupied by down and up spins, respectively 6 . Our ansatz asserts that 
any amplitude in (189) is given by 

f( Xl , ...,x n ) = E Xl ~ 1 AE X2 ~ Xl ~ 1 A ■ ■ ■ AE""- 1 ' 1 AE L ~ Xn . (191) 

Actually E and A are abstract operators with an associative product. A well defined eigen- 
function (189) is obtained, apart from a normalization, if all the amplitudes are related 
uniquely. In order to obtain the solutions through the ansatz (191) let us consider initially 
the cases of small values of n. 

n = 1. For one up spin the eigenvalue equation (188) give us three types of relations 
depending if the corresponding configuration has the up spin at x — 2, . . . , L — 1 or at the 
boundaries x — 1 and x = L, namely, 

eiE x ~ l AE L - x = -E x ~ x {E~ l AE + EAE' 1 ) E L ~ X 
- \ \{L - 5) A + h L + h R ] E x - l AE L ~ x , x = 2,...,L-l (192) 

e x AE L - x = —EAE L ~ 2 
e x E L - x A = —E L ~ 2 AE 

The solution of all these equations is obtained by identifying the matrix A as composed by 
two other matrices Bk, Ck depending on a single spectral parameter k, i. e., 

A = (B k - C k )E 2 ~ s , (195) 

with the following commutation relation with the matrix E 

EB k = e ik B k E, EC k = e~ ik C k E. (196) 



L-l 



- [(L - 3) A - h L + h R \ AE 



[(L - 3) A + h L - h R \ E A 



(193) 
(194) 



6 Differently from the solution on the periodic lattice (see (4) in section 2) in the present case it is 
not necessary to define the matrices A with the superscript s. 
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Substituting (195) in (192) and using (196) we obtain the energy in terms of the spectral 
parameter k, 

ei = -2 cos(k) - i [(L - 5) A + h L + h R \ . (197) 

Inserting (197), (195) in (193), (194) and using (196) we obtain the following algebraic 
relations 

a(k)e~ ik B k - a{-k)e ik C k = (198) 

f3(-k)e ik B k - $(k)e- lk C k = 0, (199) 
where a(k) and (5{k) are given by 

a(k) = 1 + (h L - A)e lk , p(k) = [l + (h R - A)e ik )e- tk(L+1) . (200) 
The compatibility of the relations (198) and (199) fix the spectral parameter 

- x (201) 

n = 2. The eigenvalue equation produces now fours types of relations depending on the 
relative location xi, X2 of the up spins. The amplitudes related to the configurations where 
Xi > 1, L > x 2 > x± + s give us 



e 2 E Xl - 1 AE X2 - Xl - 1 AE L - X2 = -E^- 1 [e- x AE + EAE' 1 ) E^'^E 

-E Xl - 1 AE Xa ~ X1 - 1 (E- X AE + EAE- 1 ) E 



L—X2 
L—X2 



-- [{L - 9)A + h L + h R \ E* 1 - 1 AE X ^- X AE L - X \ (202) 
2 

while the configurations where the particles are at the matching conditions X\ > 1, L > 
x 2 = Xi + 1 produce the relation 

e 2 E Xl ~ 1 AE s - 1 AE L - Xl - s = -E Xl - 2 AE s AE L ~ Xl - s 
- E Xl ~ l AE S AE L ~ x ~ s ~ l -\[{L- 5)A + h L + h R ] E Xl ~ 1 AE s ~ 1 AE L ~ x ~ s . (203) 
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Lastly the amplitudes where one of the particles are at the end points give us the following 
relations. For x\ — 1, L > x 2 > X\ + s, 

e 2 AE X2 - 2 AE L - X2 = -EAE X2 ~ 3 AE L - X2 - AE X2 - 3 AE L - X2+1 
- AE X2 - X AE L ~ X2 ~ l -\[{L- 7) A -h L + h R ] AE X2 ~ 2 AE L ~ X ' 2 , (204) 

and for 1 < x 1 < L — s, x 2 = L 

e 2 E Xl - 1 AE L ~ Xl - 1 A = -E Xl ~ 2 AE L - Xl A - E Xl AE L - Xl - 2 A 
- E Xl ~ l AE L ~ Xl ~ 2 AE -\[{L- 7) A + h L - h R ] E Xl ~ x AE L ~ Xl ~ l A. (205) 

The relations coming from the amplitudes where the particles are located at {x\ — 1, x 2 — 
1 + s) or (xi = L — s, x 2 = L) are satisfied by the solutions of (202)-(205). The solution 
of (202)-(205) is obtained by a generalization of (195), where we identify the matrix A as 
composed by n = 2 pairs of spectral parameter matrices {B kj , C^}, 

n 

A = H - C kj ) E 2 ~ s , (206) 
i=i 

obeying the following commutation relations 



EB ki = e ik *B ki E, EC ki = e^C^E 



B 2 = C 2 = B k] C kj = C k] B k] =0 j — 1, . . . ,n. (207) 



Inserting (206) into (202) we obtain the energy e 2 in terms of the unknown spectral param- 
eters hi, k 2 ] 

n y 

e n = -2j2 cos(kj) — - [(L — 1 — An) A + h L + h R ] , (208) 

where n = 2. The "bulk" relations (203) give us, after using (206) and (207), the following 
algebraic relations among the matrices {B kj ,C k] }: 

a(k h kj)B k .B kl + a(k j: k)B h B kj = 
a(-k h -kj)C kj C kl + a(-kj, -ki)C kl C kj = 

a{-k l ,k j )B k .C kl + a(k j ,-k l )C kl B k .^0 j,/ = l,...,n (j^l) (209) 
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with n = 2 and 

a(k, k') = 1 - 2Ae~ lfc ' + e -^+ fe '). (210) 

Using (206) and (207) in (204) and (205) give us the following additional relations 

a(k i )e- ik 'B k .B hl - a(-kj)e^C k .B kl = (211) 

a{-k 3 )e^C k] C kl - a(kj)e- ik 'B k .C kl = (212) 

P(-k l )e- ik "B k .B kl - (3{hV hs B k] C kl = (213) 

P{k l )e*<'C k .C*-P{-k l )e-*<'C kj B kl =0, j.l = l,...,n (j^l) (214) 

with n = 2 and 

a{k) = l + {h L -A)e ik , (3{k) = [l + (h R -A)e ik ]e'^ L+1 \ (215) 

The up to now free spectral parameters k\ and k 2 are going to be fixed by imposing the 
compatibility of the algebraic relations (209), ((211)-(214)). Using successively (211), (209), 
(213) and (209) we obtain 

ajkMkj) 2ik . {s _ 1} _ B(-k v k) 

a(-kM-k 3 f ~ B(kj,k{) ' J - X ' 2 ' ^ 3 (216) 

where 

B(k, k') = a(k, k')a(k r , -k). (217) 

General n. The eigenvalue equation applied to the components corresponding to the 
configurations where there exist no collisions (x i+1 > Xi + s; i = 1, ... ,n — 1) produces 
a generalization of the relation (202) that is solved by identifying, as in (206), the matrix 
A as composed by n pairs of spectral parameter dependent matrices {B k ,C k } satisfying 
the algebraic relations (207). In terms of the spectral parameters k 1: . . . , k n the energy is 
given by (208). The eigenvalue equation (188) when applied to the components related 
to the configurations where two particles at Xi and x i+ i are at the "colliding" positions 
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L > Xi+i = Xi + s > s give us the relations (209) for j ^ I, j, I — 1, . . . , n. The configurations 
where we have a particle at the end points x\ = 1 or x n = L give us the additional relations 



B ki B ki2 a(k h )e- ik n - C ki B ki a{-k h )e ik ^ 
C ki C ki2 a(-k h )e ik n - B ki C ki2 a(k h )e- ik n 



X k - • ■ • Xh . 



o, 
o, 



Xu - • • • X k - 



o, 

x k . ■ ■ ■ x k 



B ki _B ki j{-k ln y k ^ n ~^-^ 



C kin _C ki J{k in )e- ik ^-'^ 



0, 



(218) 



where ii,...,i n is an arbitrary permutation of the integers 1,2, ... , n, and X kj denote a 
matrix B kj or C kj . It is interesting to observe that while the algebraic relations (209) only 
relate two product of two matrix the relations (218) relate the product of n matrices. 

The matrix product ansatz (191) works only if all the amplitudes of the eigenfunction 
(189) are uniquely related. In fact the algebraic relations (207), (209) and (218) enable us 
to show that any amplitude given by the ansatz (191) is proportional to the matrix product 
B kl ■ ■ ■ B kn E L ~ n . The spectral parameters hi, . . . , k n are fixed by imposing the compatibility 
of the algebraic relations (209) and (218). For any j — 1, . . . , n we have 



B 



ki 



Bh ■ ■ ■ B k 



a(kj,ki) 



U=j+i 

,-2ikj(n-l)(l-s) 

a(k J )f3(k J 



a(k h kj) 



(3{-k j 



_ e -2ifc,-[(n-l)(l- S )-l] B 



-q a(kj,ki)a(ki, B ^ 



<*(-kj)P(-kj) /= i |^ a(k h kj)a{-kj, h 



fei 



Bh 



i ^kj -|_ i 

• • Bh 



Bh Ch 

''■ri rfj 

(219) 



that give us 



a(kj)/3(kj 



2ikj(n-l)(s+l) 



= n 



B(-kj, k 



-, j = l,...,n. 



(220) 



a(-kj)P(-kj) B(kj,ki) 

The eigenvalues are given by (208) with the spectral parameters obtained by the solutions 
of (220). The equations (208) and (220) for s = 1 coincides with the known result derived 
in [39] through the coordinate Bethe ansatz. 
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VII. CONCLUSIONS AND GENERALIZATIONS 



Formulations of a matrix product ansatz were introduced along the years in order to 
describe the ground state wave function of some special quantum chains. In general only the 
ground state wave function of these models is described by this ansatz and the quantum chain 
is not exact integrable. An exception happens in the formulation named dynamical matrix 
product ansatz where the matrices defining the ansatz are time dependent. In this case 
the full eigenspectra of some exact integrable quantum chains related to stochastic model 
are derived. We have shown in this paper that a huge family of exact integrable quantum 
chains, normally solved through the Bethe ansatz can also be solved by an appropriate 
matrix product ansatz. In our formulation, independently if the quantum chain is related 
or not to a stochastic model, the matrices are time independent. Differently from the Bethe 
ansatz where the amplitudes of eigenfunctions are given by combinations of plane waves, in 
the matrix product ansatz these amplitudes are given by a product of matrices. 

A necessary condition for the integrability of the model through the matrix product 
ansatz is the existence of an associative product among the matrices defining the ansatz. In 
addition, since the amplitudes of the eigenfunctions are related to these matrix products, the 
algebraic properties of the matrices should provide a single relation among any two matrix 
products appearing in the ansatz. These algebraic relations are obtained by imposing the 
eigenvalue equation and depend on the coupling constants defining the quantum chain. 

We have shown the formulation of the matrix product ansatz for two classes of models. 
Models with a single global conservation law (U(l) symmetry) like the XXZ chain (see 
sections 2 and 3), the spin-1 Fateev-Zamolodchikov model, and models with two conservation 
laws (U(l)<g)U(l) symmetry) like the spin-1 Perk-Schultz models, the Hubbard models as 
well the other models presented in sections 4 and 5. 

Let us discuss initially the models with a single conservation law. The associativity of 
the algebra, in this case, is immediate since the structure constants defining the algebraic 
relations among the matrices in the ansatz are complex constant numbers. In the case of 
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the XXZ chain the algebraic relations merged in the sector with two particles are enough to 
ensure the exact integrability of the quantum chain (see section 2). For the spin-1 models 
there appear additional relations involving the product of 3 and 4 matrices since the coupling 
constants defining the Hamiltonian connect up to four particles at nearest neighbor sites. 
The generalization of the matrix product ansatz for higher spin models (s > 1), although 
we did not consider in this paper follows straightforwardly. For example for spin s = 3/2 
models we should relate to the spins —3/2, —1/2, 1/2 and 3/2, in -basis, the matrices E, 
A, BE~ X B and CE~ X CE~ X C ', respectively. It will appear in this case algebraic relations 
involving the product of 3, 4, 5 and 6 matrices. 

In the case of models with two global conservation laws (sections 4 and 5) the structure 
constants defining the algebraic relations among the product of two matrices are also ma- 
trices (^-matrix) and the associativity condition is equivalent to the famous Yang-Baxter 
relations for the ^-matrices. In the case of the spin-1 Perk-Schultz model (see section 4) the 
associativity of the algebra is enough to ensure the exact integrability of the quantum chain, 
however in the case of the Hubbard model, as well the other quantum chains presented in 
section 5, there appear additional relations among the product of three and four matrices. 

Generalizations of the present matrix product ansatz to the cases where we have three 
or more conservation laws follows straightforwardly. As we have shown along this paper (see 
sections 2, 3, 4 and 6) our formulation of the matrix product ansatz allows the extension of 
several exact integrable models by including arbitrary hard-core effects, without destroying 
their exact integrability. Also it is interesting to mention that in the cases of exact integrable 
Hamiltonians associated to stochastic models, as in [17,16], since we can write all eigenfunc- 
tions in a matrix product formulation, our results imply that we can equivalently write at 
any time the probability distribution of the model in terms of a time-dependent matrix 
product ansatz as happens in the formulation of the dynamical matrix product ansatz. 

Except in section 6, all the quantum chains considered in this paper are defined on a 
periodic lattice and the eigenenergies are fixed by the cyclic property of the trace of the 
product of matrices appearing in the ansatz. In section 6 we show how to formulate the 
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matrix product ansatz in the case where the quantum chains are defined on open lattices. 
We derived the solution of the XXZ chain with magnetic fields at the endpoints of the lattice. 
The formulation of a matrix product ansatz for the other models of sections 3, 4 and 5 with 
open boundary conditions that preserve their global symmetry is also possible. 

Since the exact integrable chains considered in this paper share the same eigenfunctions 
with a related two dimensional vertex model the matrix product ansatz we formulated 
provide also a solution for these classical models. A quite interesting problem for the future 
concerns the formulation of the matrix product ansatz for the quantum chains with no global 
conservation law like the XYZ model, the 8 vertex model or the case where the quantum 
chains are defined on open lattices with non-diagonal boundary fields. 

In conclusion our results induce us to conjecture that a matrix product anstaz, along the 
lines presented in this paper, can be formulated for any exact integrable quantum chain. The 
importance of this ansatz for the future, as shown in this paper, remains on its simplicity, 
allowing quite simple generalizations and the formulation of knew exact integrable models. 
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